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SUMMARY 


A two-dimensional finite-difference numerical technique 
is presented to determine the temperature distribution of an 
internal cooled blade of radial turbine guide vanes. 

A simple convection cooling is assumed inside the guide 
vane blade. Such cooling has relatively small cooling effective- 
ness at the leading edge and at the trailing edge. Heat 
transfer augmentation in these critical areas may be achieved 
by using impingement jets and film cooling. 

A computer program is written in Fortran IV for IBM 370/165 
computer. 


INTRODUCTION 


In a turbomachine , performance benefits could be 
achieved with high turbine inlet temperature and high pressure 
ratio. Efficient expansion from a high pressure requires more 
than one axial turbine, of v/hich the first one may become 
exceedingly small. Large tip clearance, thick trailing edge 
and secondary losses in such small turbines tend to offset 
the advantages of high inlet temperatures. On the other hand, 
the small radial turbine has demonstrated an excellent per- 
formance at high stage loadings. The power output of a gas 
turbine is directly proportional to the inlet temperature 
of the combustion gases. Certainly, higher inlet temperatures 
are desired in order to achieve higher output. In recognition 
of this, researchers have devoted considerable efforts to the 
development of higher temperature blade materials and 
efficient methods of cooling. 

Radial turbine rotor cooling is accomplished by using 
cooling air passages through each rotor blade. The cooling air 
is then exhausted at the blade suction surface. The guide vane 
blades are internally cooled by streams of cooling air which 
after having cooled the vanes is ejected into the primary 
gas stream. In general, the different cooling techniques used 
for axial flow turbines, could be applied as well for radial 
turbines. For modest increases in turbine inlet temperature, 
simple convection cooled blades may be used. As further 
increases in the temperatures are required, the engine designer 
must use impingement cooling and film cooling techniques. Very 
high cycle temperature requires consideration of transpiration 
cooled blades. 

Oxidation and thermal creep are forms of failure in turbines 
with high inlet temperatures. The guide vanes are more exposed 
to failure by local oxidation than creep and hence are highly 
influenced by local hot spots. Such areas are usually localized 
at the leading and trailing edges. A theoretical method for 
estimating the temperature distribution inside the cooled guide 
vane blade, using simple convection cooling effectiveness will 
be presented. 
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GUIDE VANE COOLING 


A radial guide vane does not usually experience a high 
spanwise temperature gradient at its inlet. In addition# the 
heat conducted to the two guide vane back plates is also 
negligible. Therefore, a two dimensional solution for the 
temperature distribution is acceptable. 

The isothermal lines in an uncooled radial guide vane blade 
are shown in Figure 1 (taken from Ref. 1) . In the vicinity of 
the leading edge, the material temperature is very close to the 
gas temperature. The blade temperature drops toward the trailing 
edge. The high temperature at the leading edge region is the 
result of the high heat transfer coefficient near the stagnation 
point. 

The convection cooling requirements can be determined by 
considering the blades as heat . exchangers and evaluating the 
heat capacity of the cooling air flow through the vanes. Such 
analysis will require the solution of the heat conduction 
equation within the material combined with the heat balance 
equation for the cooling air. 

Obtaining solutions for the heat conduction equation pose 
some difficulties due to the irregular shape of the vane blade 
boundary, and the variations in the temperature and the heat 
trcinsfer coefficient around the surface. Due to the above 
circumstances, it is necessary to resort to a numerical 
solution. A finite difference technique similar to Reference 1 

t 

will be used. Modification to account for the cooling passage 
has been introduced. 

Convecrive heat transfer coefficient on both the inner and 
the outer sides of the blade surface can be estimated 
theoretically as follows: 

Hot gas side heat transfer coefficient (outside the blade) ; 

'’’he heat transfer coefficient could be obtained from the 
boundary layer characteristics around the guide vane blade. The 
boundary layer solution is obtained by using the computer program 
written by Herring and Mellor (Ref. 2) . The velocity distri- 
bution needed for the boundary layer calculations is determined 
using the computer program given by Katsanis (Ref. 3) , assuming 
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that the amount of heat transfered to the coolant from the hot 
gases will not affect the flow behavior around the blade. The 
heat transfer coefficient is then determined using the Reynolds 
analogy. Figures 2 , 3 and 4 show respectively the gas velocity, 
temperature and heat transfer coefficient distribution on the 
outside surface of the guide vane blade. 

Cooling side heat transfer coefficient (inside the vane) : 

The cooling air passes usually through axial hole passages 
of constant area, large enough to ensure turbulent flow. Thus 
the heat transfer coefficient for simple convection cooling could 
be estimated from the conventional formula given in Reference 4, 

0 8 

Nu = 0.02 Re^*“ (^) 


where 


Nu 



Re 


pVD 

U 


Tg is the material temperature on the cooling side surface 

D is the hydraulic diameter of the cooling passage. 

For other patterns of cooling the heat transfer correlation 
for flow between plates could be used, which is given as: 

Nu = 0.014 Re®*®^ 

If impingement cooling is used, the semi empirical correlation, 
given in Reference 5 could be used to determine the convective 
heat transfer coefficient. 


GOVERNING EQUATIONS AND FINITE DIFFERENCES 

In the rectangular coordinate system, the differential 
equation governing the steady state two dimensional temperature 
distriLation is given by 

i (If + - (V =0 (1) 

ax ax^ ^ 3Y aY^ ^ 

In the above equation it is assxomed that there is no heat 
generation in the blade body. Further simplifications could 


result by assxoming constant thermal conductivity, corresponding 
to the vane blade average temperature. Consequently, the heat 
equation could be reduced to: 



( 2 ) 


The heat balance equation written at a point on the outside 
surface is: 


k 


3T 

3n 


9 9 



(3) 


where n is the normal to the blade surface at the point under 
consideration . 

The heat balance equation at the blade inner surface is 
given by 




Ts) 


C4) 


Equation (2) is an elliptic equation, the boundary conditions are 
given by Equations (3) and (4) . Each boundary condition relates 
the temperature gradient normal to the surface to the local 
surface temperature. Such boundary conditions are of the 
Fourier type, relating the local heat transfer to the temperature 
on the boundary. 

To solve equation (2) , one has to use either finite difference 
or finite element technique. In this study the first approach 
is used to obtain the solution. For any grid point P, as shown 
in Figure 5, the partial differential equation C2), can be written 
in the finite difference form in terms of the temperatures at the 
points A, B, C, D and P. The points A and C are located at 
e^, £2 fractions of DX, while the points B and D are located at 
6^ and 5 2 fractions of DY from the mesh point P. For equally 
spaced grid points, the factors e^, £ 2 / 5]^ and 52 will be equal 
to 1.0. The finite difference equation is given by: 
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DX^(e, + £-) 


[— izizJ. - T. . + ■^) + — 


i/D e 


1 ^2 


- T. t|- + = 0 (5) 


DY'Uj^ + Sj) “I *'■’ "1 “2 “2 

which could be written in a more convenient form as: 



1 T. , . T. . ^ T 

i f . X/j+1 

E (£^^+£2) £2 (£^^+£2^ 


(DX) 






n 


( 6 ) 


where 


E 


_1 

1 ^2 


+ 


(DX) 

'dy^ 



~] 

2 


(7) 


The above equations are second order accurate when = 5 2 and 
“ ^2* 

At the blade boundary points, such as the point D shown in 
Figure 6, the convective boundary condition given by Equation (3) 
or (4) could be written in the following finite difference form: 


T - T h 

il ^ ^ ^ "g,c . 

an DN k g,c 


( 8 ) 


i.e. 


[T, 


= 


N 


K n DN 
+ -2jl® 


T ] 
g/c^ 


[1 + DN] 


(9) 


where 

Tjj is the material temperature at the inner or outer 
surface point D; 

Tjj is the material temperature at the interior point N; 
and g,c are subscripts for gas or cooling flow. 

Solution of Equation (6) along with Equation (9) using 
suitable iteration technique would give the temperature 
distribution in the internally cooled blade. 








COMPUTER PROGRAM 


The details of the input to the computer program are given 
in Appendix A. Initially/ the temperature matrix T is set to 
any conveniently chosen temperature. The iteration then proceeds 
using the Gauss-Seidel method. Temperatures at the interior 
points are calculated by sxibroutine "MESH" , while boundary point 
temperatures are calculated in the main progrcun. Referring to 
Figure 7, the calculation of the boundary temperatures proceed 
for each of the I lines and then in a similar manner for each 
of the J lines. The iteration process is repeated until the sum 
of the squares difference between two successive iterations is 
less than a prescribed value. 


i.e. 


l . 

all points 


N 




( 10 ) 


Where 


and 



is the temperature after the N iteration; 
is the temperature after the (N-1) iteration; 
is the prescribed error limit. 


The process of convergence could be accelerated by using 
an over-relaxation factor, a, in the following formula; 


N 


0) T^’ + (1 - 


( 11 ) 


Where 



u 


is the temperature after (N-1) iterations; 
is the temperature computed after N iterations; 
is the over-relaxation factor; 


T^ is the temperature used for the (N+l) iteration. 


The optimum over-relaxation factor is determined from different 
trials. Figure 8 gives the variation in the number of iterations, 
needed for certain convergence limit, with the relaxation factor. 

The number of iterations is significantly reduced from 440 at 
oj = 1.0 to about 140 for w = 1.6. 

After the final iteration, all the temperatures at the grid 
points are printed and the isothermals are found by linear 
interpolation . 
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The computer program could be used for simple internal 
cooling. It can also be used for internal cooling with jet 
impingement. For film cooling, the program can be used after 
some modifications to include the dimensions and configuration 
of the injection slots. 

RESULTS AND DISCUSSION 

The turbine blade and the grid network used in these 
calculations are shown in Figure 7. The blade inner and 
outer boundary surface coordinates and its material properties 
are given in Table I. 

The velocity, temperature and heat transfer coefficient 
distribution on the outer surface are shown in Figures 2, 3 
and 4 respectively. The heat transfer coefficient on the inner 
surface was determined for cooling air to hot gas mass flow 
ratio of two percent. These different data are used as input 
to the computer program to calculate the blade material 
temperature distribution. The program is listed in Appendix C 
cuid is written in FORTRAN IV for an IBM 370/165 computer. 

A sample output of the program is given in Appendix B. 

The coordinates of the boundary points through which a least 
square parabolic curve is fitted are printed in the output. 

This serves as a check on the program input parameters. The 
mean square average error is also printed for every iteration, 
making it possible to monitor the convergence process. After 
the final iteration, the internal grid temperatures and the 
surface temperatures are printed as well as the isothermal 
lines coordinates. 

The isothermal lines for the guide vane blade with single 
central cooling passages are shown in Figure 9. The difference 
between the maximum and minimum blade temperature was around 
400°K. The thermal performance of the cooled guide vane 
blades can be described by evalua+-'ng the surface local cooling 
effectiveness. The cooling effectiveness on a surface point 
is defined as: 



n 


Figure 10 shows the cooling effectiveness distribution 
along the blade surface length, for the convex and concave 
surfaces of the guid vane blade. It is clear that simple 
convection cooling has relatively small cooling effectiveness 
at the leading and trailing edges. In these areas other cooling 
methods have to be used if high local metal temperatures are 
to be avoided. At the leading edge, the cooling can be 
enhanced by application of jet impingement to provide high 
xocal cooling heat transfer coefficient. A cooling improvement 
may be achieved at the trailing edge by using film cooling. 

It can thus be concluded that simple convection cooling is only 
a primary means of cooling radial guide vanes, it must be 
enhanced by film cooling slots and impingement jets in 
critical areas. 
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D 

^ / H / C f D 
DN 

DX 

DY 

E 

h 




Hydraulic diameter 

Points indicated in Figure 5 

Distance along normal to surface of blade (Figure 3) 
(inches) 

Grid spacing in x. direction (inches) 

Grid spacing in y direction (inches) 

Constant defined by Equation (7) 

B'hu 

Heat transfer coefficient ( = ) 

hr in «F 

Grid line numbers in x direction 
Grid line n^lmbers in y direction 

^3 tU 

Thermal conductivity of material of the blade ~ o p^ 

Normal direction to the blade surface 
Nusselt Number 

A point of grid structure as in Figure 5 

Reynolds nximber of the flow inside cooling passage. 

Temperature 

Velocity inside the cooling passage 

X-coordinate 

Y-coordinate 

Density of cooling air 

Fraction of grid spacing for point B in Figure 5 

Fraction of grid spacing for point D in Figure 5 

Fraction of grid spacing for point A in Figure 5 

Fraction of grid spacing for point C in Figure 5 

Error term as defined in Equation (10) 
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Subscripts ; 


C Cooling air 

D Corresponds to boundary point in Figure 6 

g Corresponds to gas 

i,j Mesh line numbers in x and y directions 

N Corresponds to the point where normal to the boundary 

at D cuts the mesh line AP (Figure 6) 

S Surface 
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TABLE I 


PARTICULARS OF RADI2\L TURBINE NOZZLE BLADE: 


Leading 

edge Radius = 

0.0815 in 

Trailing edge Radius = 

0.0157 in 


c = 

1.94 in 

Thermal 

conductivity = 

12 

hr ft “F 

Boxindary points: 


X in 

Y (Lower) in 

Y (Upper) 

0.0 

0.100 

0.100 

0.1 

0.028 

0.130 

0.2 

0.025 

0.203 

0.4 

0.068 

0.233 

0.6 

0.105 

0.249 

0.8 

0.133 

0.255 

1.0 

0.148 

0.250 

1.2 

0.155 

0.238 

1.4 

0.147 

0.215 

1.6 

0.127 

0.183 

1.8 

0.100 

0.143 

1.94 

0.100 

0.100 

IX = O.050 in j 

DY = 0.025 in 
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FIGURE 2 


VELOCITY 


DIST. AROUND THE BLADE 




FIGURE 3. STATIC GAS TEMPERATURE DISTRIBUTION AROUND THE BLADE. 





ROUND THE BLADE 



FIGURE 5 GENERAL GRID POINT WITH 
UNEQUAL SPACrNGS 



FIGURE 6 POINTS ON A BOUNDARY 





FIGURE 7 BLADE GRID NETWORK 
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Temperatures in Degrees Rankine 
T . = Initial Temperature of the Coolant = 1400 


T^„ = Temperature of the Coolant at the Vane Cross-section » 1400 
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This computer program is an extension of the solid blade | 

program of Reference 1. The blade is drawn in the grid network | 

such that there are no multiple intersection points, i.e. each j 

grid line has no more than two intersections with the outer \ 

'i 

boundary. The blade also has to be drawn, such that 1=1 i 

i 

line is tcuigential to the blade. 5 

In order to read the coordinates of the boundary points | * 

at the grid line intersections, they are numbered in a counter- j 

clockwise direction starting from any arbiterary point of | 

intersection on the boundary as shown in Figure 7, In this | 

figure the numbering is shown only for the outer boundary; a ? 

similar procedurt can be followed for the inner boundary points. 

In case a vertical and a horizontal grid line intersect with 
the boundary at the same point, that point is given only one 
number. 

Input for the computer program; j 

The input to the program consists of four main parts; T 

1, Physical parameters for the grid system and the blade. 

2, Blade outer and inner boundary points. \ 

3, Blade surface gas temperature and cooling air temperature. 

4, Heat transfer coefficients on the outer and the inner 

boundary surfaces. 

The first part of the input consists of the following items; 

DX, DY, XK, OME, SUMM, NX, NY, Nl, N2, N3, N4, IMAX, NTE , NXO, NXI. 

These items are specified according to the format 5F5.3, 6l2, 

413, and are explained below; 

The grid spacing in the x-direction (inches) . 

The grid spacing in the y-direction (inches) . 

Thermal conductivity of the blade material (Btu/hr . in. ®F) . 

Relaxation factor (If not known, use the value 1.0) . 

The total number of mesh lines in x-direction (^ 40) . 

The total number of mesh lines in y-direction (^ 25) . , 

I mesh lines enclosing the blade inner hole without 
intersecting it; see Figure 7. In case of solid blade 
put Nl = N2 = NX. 
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DX; 

DY; 

XK; 

OME; 

NX; 

NY; 

Nl, N2 


I 


N3, N4: J mesh lines enclosing the blade inner hole without 

•intersecting it. In case of solid blade put 
N3 = N4 = NY. 


IMAX; Maximum n\mber of iterations allowed. 


NTE: Number of iterations after which a print out of the 

results is required. 

NXO: Total number of mesh lines intersections with the 

outer boundary. 120) . 

NXI; Total number of mesh lines intersections with the 

inner boundary. 120) . 


The second part of the input read in the blade inner and 
outer boundary coordinates. The Y1{I) array, corresponding 
to Y coordinates of the intersections of the I mesh lines with 
the lower surface of uuter boundary are read according to 
12F6.3 FORMAT. The next set is IP (I) point numbers (which are 
■the sequence numbers for Yl(I) points), these are read according 
to 2014 FORMAT. Similar values are read for upper surface of 
outer boundary, namely Y2(I) and IP(I). Then the Xl(J) arrays, 
corresponding to -the X coordinates of the points of intersection 
of the outer J mesh lines with -the outer boundary (12F6.3) . 

The next set is IP(J) numbers corresponding to these points. 

The arrays X2(J), IP(J) correspond to the second intersections 
of J lines with the outer boundary. All the previous coordinate 
sets are for the outer boundary. The inner boundary coordinates 
are defined in a similar manner as follows: 


Y3(I) , IP(I) 
Y4(I) , IP(I) 
X3(J) , IP(J) 


Lower surface of inner boundary 
Upper surface of inner boundary 

Near surface of inner boundary (first intersection) 


X4(J), IP(J) Farther surface of inner boundary (second intersection) 

In case of a solid blade N1 = N2 and N3 = N4, this last input 
pait should be omitted. 
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The third part of the input specifies the surface gas 
temperatures (®R) in the following order: 


TG(I,1) : 


TG(I,2) ; 


TGX(1, J) : 


TGX(2, J) : 


Gas temperature at the intersection points of the 
I mesh line with the lower surface of the outer 
blade boundary. 

Gas temperature at the intersection points of the 
I mesh lines with the upper surface of the outer 
blade boundary. 

Gas temperature at the first intersection points 
of the J mesh lines with the surface of the blade 
outer boundary. 

Gas temperature at the second intersection points 
of the J mesh lines with the surface of the blade 
outer boundary. 


The gas temperature distribution on the inner boundary are 
defined in a similar way: 

TG(I,3) Lower surface of inner boundary. 

TG(I,4) Upper surface of inner boundary. 

TGX(3,J) Near surface of inner boundary. 

TGX(4,J) Far surface of inner boundary. 


The fourth part of the input specifies the heat transfer 

2 

coefficients (Btu/hr.in °F) in a similar order as the 
temperature input: 

H(I/1) Lower surface of outer boundary. 

H(I,2) Upper surface of outer boundary. 

HX(1,J) Near surface of outer boundary. 

HX(2,J) Far surface of outer boundary. 

H(I,3) Lower ourface of inner boundary. 

H(I,4) Upper surface of inner boundary. 

HX(3,J) Near surface of inner boundary. 

HX(4,J) Far surface of inner boundary. 
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INPUT ftAT* 




^:UTM KOUffftAJIY 


IMlFH BCUNOUV 


PT.NO* . X t 


1 

0.0 

q.375000 

2 

0.007000 

0.500000 

3 

0«0260tl0 

0.673000 


O.O'dlllUll 

0.668000 

> 

Q« 062 008 

0.6SQOO0 

6 

q.lOnoilO 

0*633080 

T 

0.126000 

0.626000 

A 

a.isoouo 

0. 620000 

9 

0*2nnnnd 

0.610000 

10 

0.790000 

0.607000 

tl 

0.3VI080 

0.605000 

12 

0.39OO80 

0.601O00 

13 

O*lc0000 

0.6000OQ 

U 

O.VIPIOOO 

0*398009 

« 19 

0.630000 

0.396000 

16 

0.90QOOO 

0.389000 

IT 

0.39O000 

0.386000 

19 

0.600000 

0.380000 

19 

0.637000 

0.375000 

20 

o.65iii}go 

0*373000 

21 

0.700000 

0*366000 

22 

o»7sgooo 

0.39I080 

23 

0. 800008 

0.353O00 

26 

0.82SOOO 

0.350099 

29 

8.880000 

0.366000 

26 

0.900000 

0.336009 

2T 

0.965 100 

0*325000 

28 

0.990000 

0.326070 

29 

U8O0OHO 

0*112000 

30 

1*069999 

0. 302000 

31 

1.096999 

0*300000 

32 

1*089999 

0,290000 

33 

1*158000 

0*278000 

36 

1*1613000 

0*275000 

3S 

1.200800 

0*266000 

36 

1.265080 

0*250990 

3T 

1.269999 

C.26A0O0 

38 

1.299999 

0,211909 

39 

1*316999 

0*225000 

<•0 

1.369999 

0.213000 

61 

1*388000 

0. 200000 

62 

1*607000 

9*195000 

63 

1.6S800Q 

0*178000 

66 

1*696000 

O.175OO0 

63 

1.699999 

0.158000 

66 

1*521999 

0.130000 

67 

1.369999 

0.139000 

68 

1*382000 

0.125000 

69 

1.59199-J 

0.1 t'iifi'i 

SO 

1*661999 

0*100000 

31 

1.639000 

Q.Q9I000 

92 

t. 708000 

0.075009 

53 

1.749999 

0*036009 

56 

1,759999 

0.050000 

33 

1.799999 

o.onr^o 

96 

1.816000 

0.925004 

37 

1.869999 

0,009009 

56 

1.870000 

0.0250AQ 

39 

1.869999 

0*065000 

60 

1.861999 

0.050000 

61 

1.799999 

o.onooo 

62 

1.796000 

9.075000 

63 

1,758080 

0.100000 

66 

1.708000 

9*125000 

63 

1*650000 

0*150909 

66 

1,199999 

9*175090 

67 

1.367999 

0.200000 

68 

1,399999 

9.223090 

69 

1.695080 

0.225000 

70 

1,650000 

9*267090 

Tl 

1.64J000 

9* 250090 

72 

l,60<3a00 

0*270000 

73 

1.388000 

0*275900 

76 

1,569999 

0.291999 

73 

1.327O00 

0.300000 

76 

1,299997 

0*311999 

77 

1,266000 

0*325009 

78 

L.2 69999 

a.IIOOflO 

79 

1.200000 

9,369000 

80 

U191999 

0* 350900 

at 

1.150000 

0*366000 

82 

1.120999 

9.375090 

83 

1 .099999 

0.3P1900 

86 

1 *869999 

0.6COOCO 

83 

l.oogOQO 

0*6tf000 

86 

0*9 76000 

0.625088 

37 

0.93COOO 

0.631089 

89 

0.900008 

o*668nno 

89 

0.992088 

0.650090 

90 

0.850000 

0.662004 

91 

0.8 10800 

0.675O00 

92 

8.800000 

0.6 7*^000 

93 

0,750000 

0.6*<300Q 

96 

fl.Tnouo 

a.5caOAO 

99 

0. 700000 

0.506080 

96 

8.65>mOQ 

0.513U90 

97 

0.6IJOOOO 

4.5259O0 

98 

0.350000 

Q.5350O0 

99 

0.50)000 

0,566090 

too 

0.67'iooa 

0.650009 

101 

0.650000 

0*555090 

t02 

0.600tK10 

0,64 5001) 

103 

0.168000 

g.s7?ooo 

106 

■.1.1258811 

0. 575000 

109 

0, moOfl’J 

0.57«009 

t06 

0.25I10AQ 

0*58 7ona 

101 

0./^J7JO 

0,561800 

108 

0. 158000 

0. 696000 

109 

0. 10-1000 

0.56BOOQ 

MO 

0.060000 

0.5il2000 

til 

0 4 016009 

4,575000 

112 

0..104000 

0,550800 


PT*MO. X V 


1 

0,091980 

0.926008 

2 

0.108900 

O.**' 84 800 

1 

0.118904 

8.6190(18 

6 

O.168Q09 

0.69a080 

5 

0.151900 

0.658089 

A 

0.708800 

4.631800 

7 

0.281)900 

8.625000 

a 

9. 269000 

0.618880 

0 

i), 308800 

8.6I0400 

10 

0.3580')0 

0.607090 

11 

4.600000 

0.645804 

12 

0.650040 

0.601400 

13 

0.500000 

0.699000 

16 

0,558000 

0.391009 

15 

0.604000 

0*396999 

19 

0.613900 

0*389000 

17 

9.664009 

0*596080 

IB 

O.TOOQQC 

9.388000 

19 

0.750000 

0.375800 

20 

9.808000 

0.511499 

21 

O.flOPOOQ 

0.368880 

22 

0*058000 

0.3Ai490 

23 

0. 90090a 

8.39)009 

26 

0.925900 

0.358448 

25 

0.958000 

9.366000 

26 

0*959900 

0*336909 

27 

0.950800 

0.325000 

ZS 

9.912030 

0.326090 

29 

0.900C0O 

9.312400 

30 

O.B59C30 

O.3Q2Q0D 

31 

0.516840 

o.moooo 

32 

0. 800 COO 

0.290000 

33 

0.7*8900 

Q.ziaooo 

36 

0,739000 

0*275040 

35 

0. 709900 

9*264990 

36 

0,659040 

0.298940 

37 

0*617900 

0«Z6a949 

38 

8,600900 

0*231900 

39 

0.558990 

0.225949 

60 

9.5OO90O 

0.213909 

61 

0.6d9080 

9.298480 

62 

0*650004 

9.145Q40 

63 

8.61^0990 

0*178080 

66 

0.3MOOO 

0*179040 

65 

o.nrooo 

9.158990 

66 

0.109900 

0*150900 

67 

0*259009 

9*139000 

68 

0*7887*88 

y. 125899 

6*i 


■I. ItfOAO 

50 

O.IOODOO 

0*100990 
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REPRODUCr'uJTY OF TilE 
OBKr>J/vL FAC.i IS POOR 




Mk 


I 



DX = 0,.05000 DY = 0,02500 0MEG4= 1.60000 TKRM,CQND.= 1.00000 

NX - 39 NY = 25 MAX.ERRGR= 0.02000 

M = 2 N2 = 21 N3 = 15 N4 = 23 

NXO = 112 NXI = 50 
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tNPUT DATA 

GAS TEMP6RATUPE OISTaiUUTtON OK THE 8LA0F XNN£A AND OUTER SURFACE 


to 

KD 


o 


OnUMDARV POINTS UN l-tlNES FOR LOWER SURFACE OF 
(START FROM l» 2 TO !• 36 I 


CUTER eOUNORV 


t969*S0 197S«50 1 9 75« t 0 I975« 70 1976«70 1977.50 1976.20 1976.70 1976.90 1979.00 1979.00 1979*00 1976*60 1976*60 1976.20 
1977.60 1976.40 1975.30 1974*30 1973.45 1972.70 1972.00 1971.50 1971*20 1971.20 1971*00 1970.70 1969.90 1966.90 1968.40 
1966.00 1967.30 1966.50 1965.60 1964.70 1963.90 1963*20 


0GUN04RY POINTS ON I LINES FOR UPPER SURFACE OF 
(START FROM |« 2 TO <« 36 I 


CUTER BOUNORV 


44 1998.50 1996*50 1998.40 1996.30 1996.20 1998.10 1997.60 1997.20 1996.60 1996.60 1996.50 1996.50 1996.30 1996* 10 1995.70 

1995.20 1994.90 1994.50 1994.00 1993.50 1992.5'> 1991.50 1990.00 1966.20 1966.50 1964.60 1962.70 1980*60 1976.50 1976*50 

1974.60 1972.90 1970.60 19&9.C0 1967.50 1965.50 1963*90 


BOUNhARV POINTS ON I LINES FOR LOWER SURFACE OF 
(START FROM |« 3 TO 1 « 20 ) 


INNER BOUNDARY 


1400*00 1400.00 1400.00 1400.00 1400.00 1400.00 1400.00 1400.00 1400.00 1400*00 1400.00 1400.00 1400.00 1400*00 1400.00 
1400*00 1400.00 1400.00 


BOUNDARY POINTS ON 1 LINES FOR UPPER SURFACE OF 
(START FROM l» J TO I » 20 > 


INKER BOUNDARY 


1400*00 1400.00 1400.00 1400.00 1400*00 1400*00 1400.00 1400.00 1400.00 1400.00 1400.00 1400.00 1400.00 1400*00 1400.00 
1400.00 1400.00 1400*00 


BOUNDARY POINTS ON J LINES FOR NEARER SURFACE OF OUTER BOUNDARY 

I START FROM Ja 2 TO Js 24 I 


1963. 70 I 64 * 60 1 9 65.60 19 66 . 50 1 9 6 7 * 60 1 966 * 30 1 968.60 1970.10 1970.90 1971*20 1971.40 1972*60 1974.50 1977*00 1976*90 
1976*10 |97{*00 1972*00 1965*60 2000.00 1999.00 1999.00 1996.50 


RUUNOARV PniNTS ON J LINES FDR FARTHER SURFACE UF CUTER BOUNDARY 
(start from J» 2 TO J« 24 > 


1963*00 1964*20 1965.80 1967.50 1969*00 1970*60 1972.90 1974*60 1976.90 1979.00 1961.50 1984.00 1966.10 1966*60 1990*80 
1992.50 I993.B0 1994.60 1995*10 1996.00 1996.50 1996.70 1996.10 


BOUNDARY POINTS UN J LINES FOR NEARER SURFACE OF INNER BCUNDARV 
(START FROM Ja 16 TO J» 22 I 


14C0.00 1400.00 I4C0.00 1400.60 1400.00 1400.00 1400.00 


UnUNOARY points on 3 LINES FOR FARTHER SURFACE OF INNER BOUNDARY 
(START FROM J= 1 6 TO J- 22 I 


1400*00 |40n*0C 1400*00 1400«00 IADO.OO IAOO.OO 1400*00 




INPUr DATA 


HEAT CUfFPtCIENT OtSTRlQUTION ON THE BLADE INNER AND OUTER SURFACE 


I 

HCXJNDARV POINTS ON I-CINES FOR LOWER SURFACE OF OUTER BOUNORV 

<START FPOP !■ 2 TO 30 I 

*• 1*3a 0*66 0.96 0«99 1*00 1.01 1.00 1*02 1.09 1*06 1*12 I*!? 1*22 1*28 1*37 

1*46 |*5A U6I 1*67 1*70 1*71 1*70 1*67 1*64 1*6| 1*60 1*61 1*69 1*70 1*74 

1. 74 1*75 1.75 1* 79 1.76 1*79 1*45 


aouNOARV POINTS ON I LINES FOR UPPER SURFACE OF OUTER BOUNORV 
I START FROM |s 2 TO |« 38 ) 

•* 2*16 1*29 0*82 0*92 0*99 I«10 1.10 1.04 0.96 0.66 0*76 0*74 0*76 0*62 0*98 

0*96 0*99 1.02 1*09 1*08 1.21 1*33 1.44 1*99 1.99 1.62 1*61 1,60 1*57 1*94 

1*49 1*44 1*40 1*34 1.28 1*21 1.13 


OOUNDARV POINTS ON 1 LINES FOR LOWER SURFACE OF INNER BOUNDARY 

(START FROM 1= 3 TO l» 20 I 

0*40 0*36 0.36 0*36 0.36 0*36 0.36 0*36 0.37 0*37 0*37 0*39 0*39 0*40 0*42 

O 0*42 0*44 0.45 


eOUNOARY POINTS ON | LIKES FOR UPPER SURFACE OF INNER BOUNDARY 

(START FROM 3 TO I > 20 I 

0*40 0.36 0.36 0.36 0*36 0*36 0.36 0*36 0*37 0*37 0*37 0*39 0*39 0*40 0*42 

0.42 0*44 0*45 


BOUNDARY POINTS UN J LINES FOR NEARER SURFACE OF OUTER BOUNDARY 
(START FROM Js 2 TO Js 24 ) 

1.79 1.76 1*79 1*75 U7S 1*74 1*70 1*64 |*6l 1*62 1*67 1*71 1*66 1*50 1*20 

1.00 0*93 1*10 1*70 2.29 2.55 2.90 2*29 


( 

I 



BOUNDARY POINTS ON J LINES FOR FARTHER SURFACE OF OUTER BOUNDARY 
(START FROM Js 2 TO J* 24 1 

1*10 1*16 1*22 1*28 1.34 1.40 1.44 1*49 1.94 1.58 7.60 1*61 1.60 1*54 t*36 

1*21 1*07 1*01 0.98 0.69 0.74 0*89 1*13 


BOUNDARY POINTS ON J LINES FOR NEARER SURFACE OF INNER BOUNDARY 
(START FROM Jm |6 TU J« 22 } 

<>•44 0*42 0*37 0*36 0* 36 0*40 0*40 


BOUNDARY POINTS ON J LINES FOR FARTHER SURFACE OF INNER OOUNDARV 
(START FROM J» lb TO Js 22 I 

0*49 0*44 0*42 0*40 C*37 0*37 0*36 


I 


APPENDIX B 
PROGRAM OUTPUT 


31 


The first part of the output serves as a check for the input 
data and the calculation procedures. The coordinates of the 
boundary points, through which a parabolic curve is fitted in 
the subroutine CURVE, are printed out together with the 
coefficients of the parabola. To distinguish between the boundary 
points on the interior and exterior blade surfaces, an integer I 
is assigned the value 0 in the first case and 1 in the second case. 
To monitor the convergence process the sm of the difference 
square between two successive iterations values at all the grid 
points is also printed. The second part of the output deals 
with the temperature distribution after the final iteration. 

The temperature at all the grid points, including the fictitious 
points which lie outside the region of calculations, are printed 
out. Next, the temperature at the interior points of the blade 
are listed after eliminating the fictitious points at which the 
temperatures are constant and equal to the initialized temperature 
value. Temperature distribution on the boundary points are 
also given in the output as well as the coordinates of the 
isothermal lines. 



CU*Vl -04JI»W1‘ 


iU 

OJ 


l•«•4N|lr• amMOABy 

|» ecu>if>4*v ^niaf nubbcb 

It till BflrNi OH THf »niMn«»v* uiiCii ana am im iuaaiiimdimq antutt on aiAiit* A«t«C« *a<' ihc cmyf 

or IHI fABaniic cuAvr v-a*p«b*c*a«i 



M 

Cl 

• 1 

Cl 

4 

0.0240II 

0.41100 

O.OMOO 

0.61100 

llfl 

O.ltMOiJ 

A.1N400 
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0.19400 
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0. 42100 
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101 
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0.40000 

0.4n000 

0. )9hC0 

11 
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ITERATION* 
I TERATION* 
ITERATION* 
t TERATION* 
I TERATION* 
I TERATION* 
I TERATION* 
I TERATION* 
I TERAT tON* 
I TEPATXON* 

iteration* 

I TERATION* 
I TERATIQN* 
ITERATION* 
I TERATION* 
I TEPAr ION* 
I TERATION* 


1 

pRRCR* 

0. 63492* 

04 

2 

ERROR* 

0« 4121 IE 

04 

3 

ERROR* 

0«43164£ 

04 

4 

ERROR* 

0. 3Q796E 

04 

5 

ERROR* 

0. 28486E 

04 

6 

ERROR* 

0. 20878E 

04 

7 

ERROR* 

0. 1661 6£ 

04 

R 

ERROR* 

0* 14249E 

04 

9 

ERRCR* 

C* 12764E 

04 

10 

ERROR* 

0* 10312E 

04 

1 1 

ERROR* 

0«90S5d£ 

03 

12 

ERRCR* 

0* 76408c 

03 

13 

ERROR* 

0. 677446 

03 

14 

ERROR* 

0.S7998E 

03 

13 

ERROR* 

0* 51642E 

03 

16 

ERROR* 

0.449726 

03 

17 

ERRCR* 

0.40293E 

03 

19 

ERROR* 

0. 354336 

03 

19 

ERROR* 

0.320056 

03 

2C 

ERRCR* 

0, 283686 

03 

21 

error* 

0.2S802E 

03 

22 

ERRCR* 

0.230276 

03 

23 

ERROR* 

0. 211066 

03 

24 

ERROR* 

0. 189116 

0 3 

25 

ERRCR* 

0. 174586 

03 

26 

ERROR* 

0. 15723E 

03 

27 

ERRCR* 

0. 146046 

03 

28 

ERRCR* 

0. 131046 

03 

29 

ERROR* 

0. 12323E 

03 

30 

ERROR- 

0. 11 1536 

03 

31 

ERROR* 

0. 104846 

03 

32 

ERROR* 

0.950476 

02 

33 

ERRCR* 

3. 897826 

02 

3* 

ERROR* 

0. 9I429E 

02 

35 

ERROR* 

0. 7730 26 

02 

36 

ERR OR* 

0. 702906 

02 

37 

ERROR* 

0. 669696 

02 

38 

ERRCR* 

0.609116 

02 

39 

ERROR- 

0. S8251E 

02 

40 

ERROR* 

0.530296 

02 

41 

ERROR* 

0. 308786 

02 

42 

ERRCR* 

0. 4633S6 

02 

43 

ERROR* 

0.446876 

02 

44 

ERROR* 

0. 407376 

02 

45 

ERROR* 

0. 393496 

02 

46 

ERRCR* 

0.358466 

92 

47 

error- 

0. 346966 

02 

48 

error* 

0.315996 

02 

49 

ERROR* 

0. J0698E 

02 

50 

ERROR* 

0. 260276 

02 

51 

ERRCR* 

0.272o66 

02 

52 

ERROR* 

0. 240446 

02 

53 

ERROR* 

0. 241716 

02 

54 

ERROR* 

0.21936c 

02 

55 

ERROR* 

0. 214796 

02 

56 

ERRORS 

0« 19517c 

02 

57 

ERRCR* 

C. I9Q916 

02 

58 

ERRCR* 

0. 173886 

02 

59 

ERRCR* 

0. 170206 

92 

60 

ERRCR* 

0. 15485c 

ez 

61 

ERROR* 

0. 152I4E 

02 

62 

ERR CR* 

0. I3880E 

92 

63 

ERRCR* 

0. 1 36336 

02 

64 

ERROR* 

0. 123956 

02 


34 




T-MATRIX«GtVES ALL INTERIOR POl NTS . I N t T | AL SETTING OF T IS 1600.0 
□ALT Those insioe blade are changed and mill be different from iooo.oo 


J S 
1 

1 

2 

3 

4 

5 

6 

7 

0 

9 

10 

1 

12 

13 

14 

15 



16 

17 

14 

19 

20 

21 

22 

23 

24 

25 








I900.00 

1 900.00 

19Q0.0D 

1900.30 

19D0.0Q 

1900.00 

1990.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900,00 

1900,00 

1900,00 



1900.00 

1 900.00 

1900.00 

1900* 30 

1900. 00 

I9CO.OO 

1918.97 

1900.00 

1900.00 

1900.00 







2 

I9C0.00 

1900.00 

19CO.OO 

1900.30 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 



I9C0.00 

t 900. 00 

19C0. 00 

1900.00 

1905. ta 

1906.74 

1904.76 

1910.70 

1913,21 

19CC,00 







3 

I9O0.Q0 

1 900.00 

1900.00 

1900.00 

19C0. 00 

1909.00 

1000*00 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900,00 

1900,00 

1900.00 



1 900.00 

1 900.00 

1900,00 

1893.21 

1893. 39 

IB93.21 

1900.00 

1699.69 

1901.44 

1900.00 







4 

1900.00 

1 900.00 

1900,00 

1900.30 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900,00 

1900,00 



1900.00 

1 900.00 

1483.12 

1462.54 

1682,90 

1900.00 

1900.00 

1862.26 

1835.50 

1900.00 







5 

19C0.Q0 

1900.00 

1900. 00 

1900.00 

1900. no 

1900. DO 

1900.00 

1900.C0 

lOOO.OO 

1900.00 

1900.00 

1900,00 

1900,00 

1900.00 

1900.00 



1900.00 

1900.00 

1070.35 

1466. 13 

1900,00 

19C0.00 

1900.00 

1670.S2 

1874.35 

19CO.OO 





* 


6 

1900.00 

1900.00 

1900.00 

1900, 00 

1900. 00 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900,00 

190C,00 

1900,00 

1900.00 



1900.00 

1900.00 

14S9.20 

1656. 1 0 

1 900. 00 

1900.00 

19CO.OO 

1462*42 

1866.42 

1900.00 







7 

1900.00 

1 900.00 

1903.00 

1900.00 

1900. 00 

1900.00 

1900.00 

1900.00 

1900.00 

1900.08 

1900,00 

1900,00 

1900,00 

1900,00 

1900.00 



1900.00 

1900.00 

1050.02 

1900.09 

1900.00 

1900.00 

1900.00 

1856.64 

1660.92 

1900.00 







B 

1900.00 

1 900.00 

1900.00 

1900,30 

1900.00 

1900,00 

1900.00 

1900.00 

1900.00 

1 900* 00 

1900.00 

1900.00 

19C0.00 

1900.00 

1900.00 



1900.00 

1900.QC 

1644.31 

1900.00 

1900.00 

1900.00 

1847.40 

1851.75 

1900.00 

1 900.00 







9 

1909.00 

1 900.00 

I9C0.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900,00 

1900,00 

1900,00 

1900.00 



1900, 00 

1844.96 

1441. 40 

1900.00 

1900. 00 

1900.00 

1842.91 

1847.20 

1900, CO 

1900.00 







10 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900. CO 

1900.00 

1900.00 

1900.00 

1900,00 

1900,00 

1900.00 

1900.00 



1900.00 

1843.59 

1639.47 

1900.00 

1900.00 

1900.00 

1838.50 

1842.72 

1900.Q0 

1900.00 







11 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

19C0.00 

1900,00 

1900.00 

1900.00 

19C0.00 

1900,00 

I9on,oo 

1900,00 

1900,00 

1900.00 



1900.00 

1 843. 1 4 

19 39.27 

1900.30 

I90C.OO 

1030.41 

1034.58 

19 00. CO 

1900.00 

1900.00 






i o 

12 

1900. 00 

1900.00 

1900.00 

190D.00 

1900.00 

1900.00 

1900.00 

19CQ.00 

1900.00 

19CO.OO 

1900.00 

1900,00 

1900,00 

1900.00 

1900.00 



1900.00 

1643.61 

1439. 64 

1903.30 

1900.00 

1826.03 

1631.80 

1900.00 

19C0.00 

1900.00 







13 

1900.00 

1 900.00 

1900.00 

1900.00 

S9Q0.00 

1900.00 

1900.00 

1900.C0 

1900,00 

19C0.00 

1900.00 

19CO.OO 

1900.00 

1900.00 

1900.00 

'tj a 


1900.00 

1444.80 

1840. 76 

1900.30 

1900.00 

1627.25 

1830.46 

1900.00 

1900.00 

19C0.00 







14 

1900.00 

1 900.00 

1900. 00 

1900.00 

1900.00 

1900.00 

1900.00 

19C0.00 

1900.00 

19C0.00 

1900.00 

1900.00 

1900.00 

1900,00 

1900.00 



1050.56 

1446. 32 

1900.00 

1900.00 

1625.63 

1829.0 7 

1900.00 

1900.00 

1900,09 

1900.00 






63 2 

Itl 

19CO.OO 

1900.00 

1900.00 

1900,00 

1900.00 

1903.00 

1900.00 

190C.00 

1900,00 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

1900.00 

si 


1B53. 14 

1 448.76 

1903.00 

1900. 90 

1 929. 81 

1833.03 

1900,00 

190O.C0 

190C.OO 

1900.00 






O H 

16 

1900.00 

1 900.0Q 

1900.00 

1903.33 

i90o«no 

1900.30 

1900.00 

1900.00 

1900,00 

19CO.OO 

1900,00 

1900. 00 

1900.00 

1900.00 

1900.00 

mm 


1656. 5B 

1851.94 

1900.00 

1632.4 1 

1835.82 

1900.00 

1900.00 

1900.00 

1900. CO 

1900.00 






•jH 

IT 

1900,00 

1 900. oc 

1900.00 

1900.30 

1900.00 

1900.09 

1900.00 

1900.00 

19C0.00 

19CO.OO 

1900, on 

1900.00 

1900,00 

1900.00 

1900.00 



1060.96 

1855.66 

1900,00 

1639.69 

1842.57 

19CO.OO 

1900.00 

1900.00 

1900,00 

190C,0C 







IB 

1900,00 

1 900.90 

1900.00 

1900.00 

1900.00 

1900.00 

1900,00 

1900.00 

19CO.OO 

19CO.OO 

1900.00 

1900.00 

1900.00 

1900,00 

1673.53 



tA66.20 

1 900. on 

1847. 97 

1450. 52 

1900.00 

1900.00 

I900.00 

I904.0C 

1900.00 

1900.00 







19 

I9CO.OO 

1 900.00 

1900.00 

1900. 03 

1 900. 00 

>900.00 

1900.03 

1900.00 

1900.00 

t 900.00 

1900.00 

1900.00 

190C.00 

1900,00 

1482*97 








TEMPERATliAf- OtSTRIDliTION INSIDE THE BLADE AND ON ITS SURFACE 


tJ 

o\ 


J« 

24 

23 

Z2 

21 

20 

19 

16 

IT 

16 
15 
14 
13 
12 
1 I 
10 
9 
6 
7 
6 
5 
4 
3 


1916.14 191 3.21 1901*44 lOttS.SO 1674.35 1666*42 1660.92 1656.50 

1910*14 1910* 70 1699*69 1662*26 1670.S2 >662.42 1856*64 1051*75 1847*20 1042*72 1640*51 

1918.97 1906.76 1900.00 1900.00 1900*00 1900*00 1900*00 1047*40 1842.91 1030.50 1634.50 1631*00 1830*46 

19 1 7 * 24 1 906 * 78 1 89 3* 21 1900*00 |900*00 1900*00 1900.00 1900,00 >900*00 1900.00 1030*41 1828*03 1837*25 

1829*07 1833*03 1834,94 

I910.70 1905*18 1893*39 1682.90 1900.00 1900.00 1900*00 1900*00 1900.00 >900*00 1900*00 1900*00 1900*00 
1825*63 1829.81 1635.82 1042.57 1043.93 

1901.02 1693.23 1882. 54 1860. 13 1056,10 1900.00 1900*00 1900.00 1900*00 1900.00 1900*00 1900*00 1900*00 
1900*00 1032.41 1039*69 1850*52 1861*31 

1008*38 1883.12 1870.35 1659.20 >850.02 1644*31 1041*40 1039*67 1039.27 1039.64 1840*76 1900.00 1900.00 
1900.00 1900.09 1847.97 1661*77 1878.11 1689.19 

1846*90 1844.98 1843*59 >843*14 1843.61 1844.80 1646*32 1846.76 1651*94 1855*66 1900*00 1900*00 1675*82 
190 6*53 1926*84 

1849*94 1850*56 1651. 14 I 056. 50 1860.96 1668.20 1877.05 1900.00 >911.04 1927.74 1941.12 I94A.01 
1869.34 1873.53 1882.97 1896.69 1914.09 I92J.46 I94I.BS 1951.74 >959.17 
>899.77 1901.39 1917.64 1931,94 1943.20 1952.38 1959*50 1964.39 1965.59 
1936*56 1945.07 I9S3. 53 1959.87 1964,44 196' 58 I96B.92 
1956. 43 1960,68 1964, 85 >967.61 1969.53 .*;0.S3 
1965.46 1965.62 1967.67 1969.49 1970.52 1971.C2 
1968.60 1969.59 >970.42 1970.61 1970.88 
1970*37 1970.43 1970.66 1970.62 *970*37 
1970.55 1970.43 1970.12 1969.67 
1970,15 1969.91 1969.41 1968.75 
1969.42 1969.19 1 966* EO 1967.73 
1966, 38 1968.26 1967,49 1966* 72 
1967.24 1966, 50 1965.81 
1966.12 1965*57 I965.C9 


2 


1965.2<< 1 964.95 1964.91 


TEHPeRATURE OISfRlOUTION ON THE BOUNDARY 


\ 


HOUNOARV POINTS ON I-LINES FOR LOWER SURFACE OF CUTER OOUNDRV 
I START FROM t> 2 TO 1* 38 » 

«• 1904*30 1893*20 1883*31 1871*49 1861*14 16^2*41 1847*35 |84S*26 1644*44 1644*77 1846.06 1847*97 1650*90 1854*34 1859*09 

1865*03 1874*76 1886*21 1901*57 1919*67 1934*69 1945*92 1954*99 1961*14 1965*67 1968*17 1969*66 1970*44 1970*57 1970*33 

1969*63 1969*14 1966*23 1967*24 1966*30 1965*43 1964*91 


OOUNOARV POINTS ON I LINES FOR UPPER SURFACE OF OUTER BOUNORY 
C START FROM |a 2 TO 38 > 

4* 1914*05 1902*76 1866.78 1876*01 1667*42 1861*56 IdSS.Sl 1649*76 1843*55 1837*56 1833*21 1830.46 1830*67 1833*51 1838.16 

1842.99 1851*53 1863*30 1078*32 1904*94 1926*84 1941*04 1951*59 I9S9.6E 1964*40 1967*63 1969*64 1970*69 1971*02 1970*84 

1970*37 1969*67 1968. 75 1967*73 1966*72 1965*73 1964*93 


eOUNQARY PCIKTS ON 1 LINES FOK LUWER SURFACE OF INNER BOUNDARY 

I START FROM |a 3 TO I > 20 1 

U) 1893.73 1882.77 1867*07 1855*68 1846.41 1841*28 1638*77 1837*69 1837*65 1838*66 1840*59 1642*59 1846*04 1850*22 1855*23 

1864*34 1875*22 1892*04 


BOUNDARY POINTS ON 1 LINES FOR UPPER SURFACE OF INNER BOUNDARY 

{START FROM |s 3 TO I = 20 I 

1899*42 1880*96 |66U*74 1859.94 1653*05 1046.88 |A4|*|9 |«. 35*57 1830*07 1626*07 1824*05 1824*73 1827*20 1832*00 1837*90 
1647*41 1860*01 1873*27 


boundary points on J lines for nearer SURFACE OF CUTER BOUNDARY 
{START FROM Jm 2 TO J- 24 | 

1965*20 1966*12 1967.24 1968*38 1969*42 1970.15 1970.55 1970*37 1968*60 1065*46 1956*43 1936*56 1899*77 1869*34 1849*94 
1846.90 1888*38 1901*02 1910*70 1917*24 1918*97 1918*14 |916*|4 


BOUNDARY POINTS UN J LINES FOR FARTHER SURFACE OF CUTER BOUNDARY 
{START FROM 3= 2 TO J= 24 1 

1964*91 1965*09 1965.8J 1966.72 1967*73 1968.75 1969*67 1970*37 1970*88 |9T1*02 1970*53 1968.92 1965*59 1959*17 1946*01 
1926*84 1889*19 1661*31 1843*93 1834*94 1R3C.46 1840*51 1858*50 


BOUNDARY POINTS ON J LINES FOR NEARER SURFACE OF INNER BCUNOAPV 
{START FROM 16 13 22 I 

1882*88 1856*28 1841*03 1049*94 1682* 58 1868*71 |900*16 


BCXiNOARY POINTS ON J LINES FOR FARTHER SURFACE OF INNER BOUNDARY 
I START FROM 16 TO Ja 22 I 

1901*56 1863*45 1843*87 1630*81 1824*05 IH3I.01 1849*23 




w'« 


ISO^HSRMAC LINS LOCATIONS 


1 

J 

T 

T - t 

T - 2 

PRAC 

2 

19 

1900 

1904, 3030 

19C5. 1 763 

C.7982 

2 

20 

1906 

1905. 1 763 

1906. 7817 

0.5131 

2 

21 

1907 

1906.7817 

190 6. 7649 

0.1101 

2 

21 

100 8 

1906.701 7 

19C3. 7649 

0.6143 

2 

22 

1909 

1908. 7649 

1910.6960 

0.1217 

2 

22 

1910 

1900.7649 

19;0.6960 

C.6396 

? 

23 

191 1 

1910. 6960 

1913.20 85 

".laiO 

2 

23 

1912 

1910.6960 

1913.2085 

O.S190 

2 

23 

1913 

1910.6960 

191 3.2035 

0.9170 

2 

2i 

1914 

1913.2006 

1914.0481 

C.9427 

3 

23 

1900 

1899.6077 

190 1.4368 

C. 1 705 

3 

23 

1901 

1899.6677 

1901.4368 

0.7503 

3 

24 

1901 

1901.4360 

1900.0000 

0.3040 

A 

10 

1003 

1883. 1 194 

1862.5437 

0.2074 

4 

22 

laai 

1880.9595 

1 362.2837 

0.0306 

4 

22 

1 892 

iaeo.9595 

1382. 2837 

0*7858 

4 

23 

1093 

1082.2637 

1835. 5039 

0.2224 

4 

23 

1 894 

1002.2837 

1305*5039 

0.5330 

4 

23 

1880 

1882. 2837 

leaS. 5039 

C.0435 

4 

24 

1 986 

1883.5039 

19JO.OOOO 

0.0342 

5 

19 

1869 

I860. 1252 

1667. 0676 

0*1 184 

S 

la 

1069 

10 70.3477 

1363. 1252 

0.6064 

3 

10 

1079 

1870.3477 

1663.1252 

0.1 564 

5 

17 

1971 

1071. 4066 

1870. 3477 

C.4272 

5 

22 

1869 

1060*7393 

1870. 5227 

0* 1462 

3 

22 

18 7C 

1860.7393 

1670. 5227 

0*7069 

m 

23 

1071 

1070.5227 

1874.3472 

0.1 243 

5 

23 

1 872 

18 70.5227 

1874.3472 

0.3363 

5 

23 

1973 

1870* 5227 

1874, 3472 

0.6477 

5 

23 

1974 

1070, 5227 

1874, 3472 

C.9092 

5 

24 

1975 

1874.3472 

1900.0000 

0.0254 

5 

24 

1875 

1874. 3472 

1 900.0000 

G.C644 

6 

19 

1 956 

1056.0977 

1655. 6829 

C.2354 

6 

Id 

1957 

1859.2043 

1856.0977 

0.7095 

6 

1 a 

1858 

1659. 2043 

1 356. 0977 

0.3077 

6 

1 a 

1 859 

1859.204J 

1356.0977 

C.0656 

6 

17 

1060 

1061. 1 382 

1359* 2043 

0.5386 

6 

17 

1061 

1361. 1 382 

1359. 2043 

0.071S 

6 

22 

I960 

1859.9350 

1302.4224 

0.0256 

6 

22 

1861 

1859, 9350 

1862. 4224 

0.4200 

6 

22 

1862 

1859.9353 

1362.4224 

0.3301 

6 

23 

1963 

1062.4224 

1366.4153 

C.1447 

6 

23 

1864 

1062. 4224 

1866. 41 53 

n,395l 

e 

23 

1065 

18 62.4224 

1866.41 53 

C.6456 

e 

23 

1966 

1862. 4224 

1866. 41 53 

0.8960 

6 

24 

1 867 

1866.4153 

1 900* 0000 

0.0174 

7 

18 

1947 

1350.0190 

1 346.4099 

C .6365 

7 

Ifl 

1840 

1850.0190 

164 6.4099 

C.5094 

7 

10 

184Q 

1050.0 190 

1846.4099 

C.2024 

7 

18 

1 850 

1 050. 0190 

1846.4099 

C.D053 

7 

17 

1051 

10S2.4 133 

1 650. 019C 

0.5904 

7 

IT 

1 852 

1052.4 133 

1850.0190 

C. 1720 

7 

22 

1054 

1053. 0510 

1356. 6367 

0.2647 

7 

2? 

1055 

1353.0510 

1356.6367 

0.5*15 

T 

22 

1 856 

1053.051 0 

1856.6367 

6.3224 

7 

23 

1 057 

18 56.6 36 7 

I860. 91 60 

0.0349 

7 

23 

1850 

1856.6 307 

1360.91 60 

0.3136 

7 

23 

1 859 

18 56.6367 

1360.9I6C 

0.5523 

7 

23 

1 960 

1056.6367 

I860. 9 1 6C 

C.7059 

7 

24 

1 861 

1860.9 160 

1900.0000 

0.0C21 

e 

la 

194 2 

1044.3061 

I84t. 2706 

0. 7619 

8 

18 

1843 

1844.3061 

ia«l .2730 

0.4310 


38 



APPENDIX C 


PROGRAM LISTING 


o o o o o o 


:v G LEVEL 21 


(«A IN 


DATE = 74346 


16/55/51 


COMMON T(40,25) *TB(40,4) »TBX(4t25) t TG (40 »4 ) , TGX ( 4 , 25 ) , H ( 40 , 4 ) t HX { 4 
lt25) ,Y1 (40) ,Y2(40) ,Y3( 40J tY 4(40) ,X1( 25) ,X2( 25) ,X3( 25) ,X4(25) ,BX( 12 
20) ,BY( 120) »BX1 ( 120) tBYl ( 120) » IP( 120 ) 

COMMON OX,OY,NX,NY ,NX0,NXI ,1 TER,N1 tN2,M3 ,N4 

REA0(5,5) OX,OYtXK,CME,SUMM,NX,NY ,N1 ,N2f N3,N4, I MAX , NTE , NXO , NX I 
5 F0RMAT(5F5.3,6I2,4I3) 

NX1=NX-1 
NY1=NY-1 
N11=N1+1 
N21=N2-1 
N31=N3+1 
N41=N4-1 
DO 46 I=1tNX 
DO 46 J=1,NY 
46 T(I,J)=1600. 

DO 44 1=1, NX 
DO 44 J=l,4 
44 TB( I,J )=1600. 

DO 49 J=1,NY 
DO 49 1=1,4 
49 TBX( I, J )=1600. 


READS IN BOUNDARY POINTS SET AND THE CORRESPONDING SEQUENCE NUMBERS 
FOR THE boundary POINTS 


READ(5,10) ( Yl( I ) , I = 1,NX1 ) 

10 FORMAT! ( 12F6.3) ) 

REA0(5,11) ( IP(I),I = 1,NX1) 

11 F0RMAT(20I4) 

DO 12 1=1, NXl 
X=FLOAT( I-1)*DX 
BX( IP( I ) )=X 

12 BY( IP{ I ) )=Y1( I ) 

REA0(5,10) ( Y2( I ) ,I = 1,NX1 ) 

READ(5,11) ( IPd ) ,I = l,NXl ) 

DO 16 1=1, NXl 
X = FLOAT( I-l )’!'DX 

BX( IP! I ) ) =X 
16 BY( IP( I ) )=Y2( I ) 

READ(5,10) (XKJ) ,J = 1,NY1) 

READ(5,11) ( IP( J ) ,J = 1,NY1) 

DO 20 J=1,NY1 
Y=FLOAT( J-1 )*DY 
BX(IP( J))=X1( J) 

20 BY! IP! J ) ) =Y 

READ(5,10) !X2! J) ,J = 1,NY1 ) 

READ(5,11) ! IP! J) ,J=1,NY1) 

DO 24 J=1,NY1 
Y=FL0AT!J-1)*DY 

BX(IP!J))=X2(J) REPROOUCiSail’Y OK 1 HK 

24 BY!IP(J))=Y QjaeiHAb PAGE IS POOR 
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IFIN1.EQ.N2.AND.N3.EC.N4) GO TO 80 
RE AO (5, 10) (Y3( I ) f I =Nll,N2l) 

READ(5,U) ( IP(I ), I = N1L,N21) 

00 28 I = M1,N21 
X==FLOAT ( I -1 ) ’I'OX 
Bxi { ip( I) : -X 
28 BYI ( IP ( I) )=Y3(I ) 

READ(5,10) I Y4( I ) , I=Nl I, N21) 

REA0(5,11) ( IPII )f I = N11,N21) 

DO 32 I=M1,N21 
X=FLOAT ( I-l) «0X 
BXI { IP( I ) ) = X 

32 BYI ( IP{ I) )=Y4( I) 

REA0(5,10) (X3IJ ), J = N31,N41) 

READ{5,11) ( IPIJ)» J=N'3ltN41) 

DO 36 J = N31,N41 
Y=FLOAT( J-l )*DY 
BXI IIP! J) )=X3( J) 

36 8YH IP( J) )= Y 

READI5,10) (X4(J), J=N31,N41) 

READ{5,11) ( IPU ), J = N31,N41) 

DO 40 J=N31,N4l 
Y = FLCAT( J-D’i'OY 
BXI ( IP( J) )=X4( J) 

40 BYI { IP( J ) )=Y 
80 CONTINUE 

READS IN GAS TE«PS AND THEN HEAT TRANSFER COEFFICIENTS 

REA0(5 ,50HTG{ 1,1 ) , 1=1 ,NX1 ) 

READ(5,50)(TG{ I,2)tI=l,NXl) 

IFIN1.EQ.N2) GO TO 82 
READ! 5, 50) ( TGI 1 , 3 ) , I =N 1 1 ,N2 1 ) 

REA0(5 ,50 J (TGII,4) , I=N11 ,N21) 

82 CONTINUE 
50 FORMAT ( 10F8.2) 

READ! 5 ,50 ) ITGXI 1, J ) , J=2,NY1) 

READI5,50) (TGX(2,J) ,J=2,NY1) 

IF(N3.EQ.N4) GO TO 33 
READ(5,50) (TGX(3, J ) ,J=N31,N41) 

READI5,50) (TGXI4, J) ,J=N31,N41) 

33 CONTINUE 

RE AD (5,45) I HI I ,1 ) ,1=1, NXl) 

READ! 5,45 ) (HI I ,2 ), 1 = 1, NX 1 ) 

IFI N1.EG.N2) GO TO 84 
REA0I5,45) (HII ,3) ,I=M1,N21) 

READ(5,45) IH(I,4), I = N11,N21) 

84 CONTINUE 
45 PQRMAT I 10F8.4) 

READI5,45) (HXI 1, J ) , J=2,NY1) 

REA0I5,45 )IHXI2, J) , J=2,NY1) 

IF (N3.Eg.N4) GO TO 85 
READ(5,45 HHXI 3,J ), J=N31,N41) 
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READ(5,45 ) (HX(4, J), J=N31»N4l J 

85 CONTINUE 
ITER=1 
WRITE(6,86) 

86 FORMAT! // ,25X, • INPUT DATA') 

WRITE(6t88 ) 

88 FORMAT!///, 15X, 'BLADE BOUNDARY POINTS COORD INATES* /////17X , 'OUTER 
1 BOUNDARY' ,9X, • INNER B CUNDAR Y ' , / / lOX , » PT . NO . • , 5X , • X • , 9X , ' Y • , 6X , 
2'PT.N0.',6X, 'X',9X, 'Y'/) 

NXI1=NXI+1 

WRITE!6,851) ! ! I, BXI I ),BY! I ) , 1,8X1 1 1 ) ,BY! I) ) , I = 1,NXI ) 

» • WRITE!6,8 53) !! I,BX!I ), BY! I)) , I=NXIl,NXO) 

j 851 F0RMAT!9X,I6,2F10.6,I6,2Fia.6) 

? 853 F0RMAT!9X,I6,2F10.6) 

\" WRITE! 6,561) 

' 561 FORMAT !• 1' ,//15X, • INPUT DATA •// 1 5X ,' GAS TEMPERATURE DISTRIBUTION 

' ION THE BLADE INNER AND OUTER SURFACE',//) 

WRITE!6,530) NXl 
WRITE!6,525) ! TG ! I , 1 ) , I = 1 ,N XI ) 

WRITE!6,531) NXl 

: WRTTE(6,525) ! TG II , 2 ) , I = 1 , NX 1 ) 

i IF!N1.E0.N2,AND.N3.E0,N4) GO TO 98 

' WRITE 16,532) N11,N21 

WRITE!6,526) 1 TG ! I , 3 ) , I =N1 1 , N21 ) 

WRITEi6,533) N11,N21 

WRITE! 6,526) I TGI 1 , 4 ) , I = N 1 1 , N21 ) 

98 CONTINUE 

; WRITE 16,534) NYl 

i WRITE!6,526) ( TGX! 1 , J > , J =2 ,N Y 1 ) 

' WRITE!6,535) NYl 

WRITE!6,526) ! TGX ! 2 , J) , J=2 ,NY I ) 

IF!N1.EQ.N2.AN0.N3.EC.N4) GO TO 96 
WRITE!6,536) N31,N41 
WRITE!6,526) !TGXI3,J) ,J=N31,N41) 

WRITE!6,537) N31,N41 

WRITE! 6,526) 1 TGX ! 4 , J) , J = N31 ,N41 ) 

WRITE! 6,401) 

96 CONTINUE 

j WRITEI6,562) 

562 FORMAT! ' 1* ,//15X, » INPUT DATA '// 1 5X ,' HEAT COEFFICIENT DISTRIBUTION 
ION THE BLADE INNER AND OUTER SUR=ACE’,//) 

WRITE!6,530) NXl 
WRITE!6,525) I H! I , 1) , 1 = 1 ,NX1 ) 
j' WRITE!6,53l) NXl 

i WRITE!6,525) ! H! I , 2 ) , I = i , NX I ) 

* IF! N1 .EQ.N2.AND.N3.EQ.N4) GO TO 97 

I WRITE 16,532) N11,N21 

WRITE! 6,526) ! H ! I , 3 ) , I =N11 , N2 1 ) 

WRITE16,533) Nll,N21 
WRITE!6,526) I H I 1 , 4 ) , I = N1 1 , N21 ) 

97 CONTINUE 

WRITE !6,534) NYl 
WRITE!6,526) I HX ! 1 , J ) , J = 2 ,NY 1 ) 
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W3ITE(6r535) NYl 
WRITE(6,526) ( HX ( 2 t J ) t J = 2 i N Y 1 ) 

I^(N1.EQ.N2.AN0.N3.EQ,N4) GO TO 95 
WRITE(6,536) N31tN41 
WRITE(6,526) ( HX ( 3 , J ) , J= N3 1 »N41 ) 

WRITE(6,537) N3l,N41 
WRITEI6»526) ( HX ( 4 , J ) , J= N3 I ,N41 ) 

WRITEI6,401) 

95 CONTINUE 

PRINT TITLE FOR SUBROUTINE CURVE OUTPUT 
WRITE(6,43) 

43 FQRMAT(1H1,//,5X, 'SUBROUTINE CURVE -OUTPUT- ' //5X 1=0, CUTER BOUN 
IDARY. 1=1, INNER B OUNDA RY • / 5X , • I P BOUNDARY PCINT NUMBER') 

WRITEI6,41} 

41 FORMAT! /,5X,' (B2,C2 ) IS THE POINT ON THE BOUNDARY. (Bl,Cl) AN 

ID (B3,C3) ARE THE SURROUNDING POINTS ON BLADE. A,B,C, ARE THE COEF 
2S'/5X ,'OF THE PARBCLIC CURVE Y= A +B*X+C* X*X • , / / , 4X , ' I • , 3X, ' I P • , 5 X , 
3'Bl',gx,'Cl',3X,'B2',8X,'C2’,8X,'B3',8X,'C3',9X,*A',9X,'B',9X,'C', 
4// ) 

1 CONTINUE 

STARTING OF GAUS S-S E I DEL L METHOD OF ITERATION ON I LINES 


SUM=0. 

on 250 1=2, NXl 

IL = IFIX(Y1 ( I )/DY+0.0001) +2 

I H=IFI X( Y2 ( I) /DY-l-0.000 1) +1 

IFI ABS(Y2( I ) -FLOAT! lH-1 )=^DY) .LT. 0.00001) IH=IH-1 

OUTER LOWER BOUNDARY 

B2 = FL0AT! I-1)«0X 
C2 = Y1 ( I ) 

CALL CUR(B2,C2,B,C,0) 

XMN=B+2. *C#B2 
IF! I.EQ.NXl ) XMN=0.01 
YT=FLOAT! IL-1)*DY 
IF! ASS(XMN) .LE.0.015)G0 TC325 
X= (C2-YT )#XMN+R2 
IF!X.LE.82)G0 TQ315 
XR=B2+DX 

IF(!X2!IL)-XR) .GT. 0.0001! GO TO 324 
TTT = T{ I, !L) + !TBX!2 , 1 L ) -T { I , I L )) -S' I X -82 )/ ( X2 ! IL ) -B2 ) 
GO TO 320 
324 CONTINUE 

TTT=T !I,IL)+(T(H-1,IL)-T!I,IU )*(X-B2)/DX 
GO TO 320 

315 TTT = T! I , I L ) -! T ( I , I L )- T ! I -1 , I L ) ) * ( B2-X ) / DX 
320 XNL=! X-32)=f'(X-B2) + (C2-YT)«!C2-YT) 

XNL=SORT( XNL ) 

GO TO 330 
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325 TTT=T(I»IL) 

XNL=YT-C2 

330 AL=H(I,n/XK 

TNEW=t XNL=^AL*TG{ I ,n+TTT)/I l.+XNL«AU 
TNEW=( l.-0ME)*TB( I , 1 ) +nHE*TNEW 
SUM=SUM+(TB(I ,l)-TNEW)«{Te( I,1)-TNEW> 
TB( 1,1 )=TNEW 

IFn.GT.Nl.AND.I.LT.N2 ) GO TO 230 


INTERIOR POINTS STARTING FROM LOWER END FOR I LINES WHICH DO NOT CUT 
THE INNER 0CUNCARY 

00 210 J=IL,IH 
CALL MESHI I ,J,TNEW) 

TNEW=( l.-OME)*T{ I , J ! +OME# TNE W 
SUM=SUM+{TNEW-Tt I , J )) * ( TNEW-T ( I,J ) ) 

210 T(I,J»=TNEW 
GO TO 249 


FOR I LINES INTERSECTING INNER BOUNDARY-ALL INTERIOR POINTS FROM 
LOWER SURFACE OF OUTER BOUNDARY TO LOWER SURFACE OF INNER BOUNDARY 

230 CONTINUE 

ILI=IFIX(Y3( n/0Y+0.0001)+l 

IFIABSI Y3( I )-FLOAT t ILI-I)*DY),LT.0.00001 ) ILI=ILI-1 

IHI=IF1X( Y4( I I/OY+O.QOOl )+2 

DO 235 J=IL,ILI 

CALL MESHI I , J.TNEWi 

TNEW=( 1 .-OME)=«‘T( I , J)+OME*TNEW 

SUM = SUM+(TNEW-T( I , J ) ) * I TNEW-T { I , J )) 

235 T(I,J)=TNEW 

INNER BOUNDARY - LOWER SURFACE 

B2=FL0AT( I-l )»DX 
C2=Y3( I ) 

CALL CUR(B2,C2,3,C, 1) 

XMN=B+2.*C*B2 
IF{I,E0.N21) XMN=0.01 
YT=FLOAT( ILI-ll’i'DY 
IF(A8StXMN),LE.0.015) GO TO 355 
X=( C2-YT) X£XMN+B2 
IF(X.LE.B2J GO TO 350 

TTT=TI I,ILI)+(T(I+1,ILI)-TII,ILI) J*IX-B2)/DX 
GO TO 352 

350 TTT=Tl I ,I LI S-I TI I, I LI)-TI I-U an }*' 62-X) /OX 
3 52 XNL = (X-B2)’!‘{X-32) + {C2-YT) = (C2-YT) 

XNL=SORT( XNL ) 

GO TO 358 
355 TTT=Tt I ,I LI) 

XNL=C2-YT 
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358 AL=H( I »3) /XK 

TNEW=IXNL*AL*TG( 1,3 J+TTT )/ U.+XNL*AL) 

TNEW=( 1,-0ME)«TB( I ,3 ) + Q^^E'!'TNEW 
SUM=SUM + (TB( I,3)-TN£W»*( TB ( I ,3)-TNEW) 
TB(I,3)=TNEW 

UPPER SURFACE OF INNER BOUNDARY 
C2=Y4( I ) 

CALL CUR(B2,C2,B,C,1) 

XMN=B+2.*C*B2 

IF(I.EQ,N21) XHN=0.01 

YT=FLOAT{ lHI-1 )*0Y 

I = UBS( XMN).LE. 0,015) GO TC 370 

X=(C2-YT)#XMN+B2 

IF(X,LE,B2) GO TO 360 

TTT=T(I ,THI)+(TI I + 1,IHI)-T (I ,IHI) )*(X-B2 )/DX 
GO TD 365 

360 TTT=T( I, IHI)-(T( I, IHI )-T( I-l, IHI ) )*(B2-X) /DX 
365 XNL=(X-B2)*(X-B2)+(C2-YT )«( C2-YT) 

XNL = SQRT( XNL ) 

GO TO 375 
370 TTT=TI I , IHI ) 

XNL=YT-C2 

375 AL=H(I,4)/XK 

TNEW=(XNL*AL*TG( I ,4 ) +T TT ) / ( 1 . +XNL*AL ) 

TNEW=I 1 ,-CME)*Ta( I,4)i-CME>«'TNEW 

SUM = SUM+(TB( 1,4)-TNEW)’^(T8(I ,4)-TNEW) 

TB(I,4)=TNEW 


INTERIOR POINTS FROM UPPER SURFACE OF INNER BCUNDARY TO THE UPPER 
SURFACE OF OUTER BOUNDARY 

DO 240 J=IHI ,IH 
CALL MESHU ,J,TNEW| 

TNEW=t l.-OME)«T( I, J) +CME*TNEW 
SUM=SUM+( TNEW-TI I , J ) ) # ( TNE W-T ( I , J ) ) 

240 T(I,J)=TNEW 
249 CONTINUE 

UPPER ROUNOARY OF CUTER SURFACE 
C2=Y2( IS 

CALL CUR(B2,C2,8,C,0) 

XMN=e+2.*C*P2 
IF(I.EO.NXl) XMM=0.01 
YT = FLOATl IH-1 )»0Y 
IF{A8S(XMN).LE.0.0.5S GO TC 340 
X= (C2-YT )*XMNfB2 
XL=B2-DX 

IF((XL-X1 (IH)).GT.O 00001) GO TO 336 
TTT=T( I, IH)-(T( I, IH)-TBX ( 1, I H ))*( B2-X )/( B2-X 1 ( IH) ) 


BEPRODUClBtUTV OE^UE 

OKBIWAt PA08 IS fOO* 
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GO TO 337 
336 CONTINUE 

IFIX.LE.B2) GO TO 335 

TTT=T( I,IH)+(T( I + 1,IH)-T( I , I H )) * ( X-B 2 ) /DX 
GO TO 337 

335 TTT=T( I , I H )- ( T ( I , I H )-T ( I- 1 , I H ) )*(B2-X)/0X 
337 XNL=(X-B2)#(X-B2)+(C2-YT J*(C2-YT) 

XNL=SQRT( XNL) 

GO TO 345 
340 TTT=TII,IH} 

XNL=C2-YT 
345 AL=H(It2)/XK 

TNEW=( XNL*AL*TG( 1,2 )+TTT )/ll.+XNL*AL) 

TNEW={ l.-OME)*TB( I , 2 ) +OME*TNEW 
SUM=SUM + {T8( I,2)-TNEW)>!‘ITB(I ,2)-TNEW) 

TB( I,2)=TNEW 
250 CONTINUE 

ITER. FOR J LINES ONLY BOUNDARY POINTS INTERSECTED BY J LINES 

DO 490 J=2,NY1 
IN =IFIX{X1U)/0X + 0.0001 ) + 2 
lA =IFIX(X2{ J)/OX+O.OOOU + 1 
IF(A9S(X2tJ )-FL0AT( IA-1)*CX) . LT.O. 00001 J IA=IA-1 

NEAR SURFACE OF OUTER BOUNDARY 

B2=X1( J) 

C2=FL0ATIJ-1)*DY 
CALL CUR (B2,C2,B,C,0) 

XMN=B+2.#C*82 
XT*FLOAT( IN-1 )*0X 
IF( ABS(XMN) .GE.5. )G0 TO 410 
XL=XT-OX 
X=-0Y*XMN+82 
IF(X.LE.B2) GO TO 406 
IFl(X2(J^-l)-XTJ.GT, 0.0001) GC TO 405 
TTT=T{IN-1 ,J+1H-!TSX(2, J+l)-T(IN-l, J+1) )*! X-XL ) / I X 2 1 J+ 1 )-XL) 
GO TO 407 

405 CONTINUE 

TTT=T{ IN,J+1)>( T( IN-1 ,J+1 )-T( IN,J+1 ))*{XT-X)/DX 
GO TO 407 

406 CONTINUE 
X=DY*XMN+B2 

TTT=T( IN,J-l) + (TBX( 1,J-1)-T( IN,J-1) > « ( XT-X ) / ( XT- X 1 { J- I ) ) 

407 XNL=IB2-X)’!'{32-X}+DY*DY 
XNL=SQPT( XNL) 

GO TO 415 
410 TTT=T(IN,J) 

XNL=XT-B2 

415 AL=HX( 1 , J ) /XK 

TNEW-(TTT + AL*XNL*TGXf ’ ,J) )/( l.+XNL*AL) 

TNEW=( l.-OME )’*'TBX( 1,J) tOME^TNEW 
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SUM=SUM+( TMEW-TBXI 1, JJ TNE W-TBX( U JJ j 
T3X( 1, J)=TNEW 

IF|J.GT,N3.AND.J,LT.N4) GO TO 418 
GO TO 445 


FOR J LINES INTERSECTING INNER 8GUNDARY-BQUN0AR Y POINTS ON THE 
NEARER SURFACE OF INNER BOUDARY 


418 CONTINUE 

INI=IF IX( X3( J)/DX+0,0001 ) +1 

IF( ABS( X3 IJ) -FLOAT ( INI-U*DX ). LT , 0.0000 1 ) IM=INI-1 

IAI=IFIX(X4( J)/DX+0,0001)+2 i 

B2=X3IJ) *] 

CALL CUR(B2,C2,B,C,1) 

XMN=B+2.=«=C«e2 
XT=FLOAT( INI-1)*DX 

IF(ABS(XMN) .GE.5. )GC TO 425 ] 

X=DY«XMN+B2 

IF(X.GT.R2) GO TO 421 | 

TTT=T( lNTfJ-l) + (T{INI+ltJ-U-TlINI,J-l) J«(X-XT)/DX j 

GO TO 422 I 

421 X=-DY*XMN+B2 I 

TTT=T( INI ,J + 1)+(T( IM, J + D-T ( INI+1,J+1) )’i'(X-XT)/DX : 

422 XNL=IB2-X)*(B2-X)+DY*DY j 

XNL=SQRT { XNL) 

GO TO 426 

425 TTT=T(INI,J) 

XNL=B2-XT I 

426 AL=HX( 3» J ) /XK 
TNEW=(TTT+AL*XNL*TGXI3tJ) )/( l.+XNL*AL) 

TNEW=( 1 .-CME)*TBX{ 3, J) +OME*TMew 

SUM = SUM + ( TNEW-T8XI 3f J) J * ( TNE W-TRX ( 3 t J ) ) 

TBX(3f J )=TNEW 

FARTHER SURFACE OF INNR BOUNDARY 
32=X4(J) 

CALL CUR(B2»C2 ,B,C,1) 

XMN=B+2.*C*32 
IF(J.EQ.N31) XMN=6.0 
XT=FLOAT( I AI-1)*DX 
I = ( ABS(XMN) .GE.5. )GC TO 435 
X=-DY’i'XMN+B2 

TTT = T{ lAI , J+1 )f (T( IAI-1, J + 1 )-T < I AI , J+1 ) )*( XT-X)/nx z' 

430 XNL=f 82-X)# {B2-X)+DY*DY j 

XNL=SQRT( XNL) J 

GO TO 440 3 

435 TTT=T(IAI,J) ^ 

XNL=XT-H2 j 

440 AL=HX (4, J )/XK ] 

TN5W={ TTT+AL*XNL*TGX(4f J) ) /( l.+XNL^AL) 3 
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TNEW=( l.-0ME)*TBX(4, J) +OME+TNEW 
SUM=SUM+(TNEW-TBX(4, J) )« ( TNE W-T8X ( 4» J) ) 

T3X(4, J)=TNEW 

FARTHER BOUNDARY POINTS ON OUTER SURFACE 

445 B2=X2IJ) 

CALL CUR(82,C2fB,C»0) 

XMN=B+2.#C*B2 
IF(J.EQ.2 ) XMN=6.0 
XT=FLOAT( IA-1)*DX 
I<=(ABS(XMN) .GE.5.) GO TO 455 
XR=XT+DX 
X=DY*XMN+B2 

IF( (XT-XK J-1) ).GT. 0.00001) GO TO 450 
TTT=TBX{ l,J-l)-(TBXa»J-l)-T( IA + 1,J-1) )*(X-X1(J-1) ) / ( X R-X 1 ( J- 1 J } 
GO TO 452 
450 CONTINUE 

TTT=T( lA , J-l)+( T( IA+1 ,J-U-T(I At J-1 ) )*( X-XT)/DX 
452 XNL=l B2-X)«(B2-X)+CY*CY 
XNL=SQRT( XNL) 

GO TO 460 
455 TTT=T(IA,J) 

XNL=B2-XT 

460 AL=HXI 2t J )/XK 

TNEW={ TTT + AL*XNL*TGX(2»J) ) / ( 1 . <-XNL«A L ) 

TNEW=( 1,-OME)*T0X( 2, J) +OME*TNcW 
SUM=SUM + {TBX(2t J)-TNEW ( T8X ( 2 t J )-TN EW ) 

TBX(2t J)=TNEW 
490 CONTINUE 

DO 495 I=ltNXl 
DO 495 J=ltNYl 
X=FLOATn-l)*OX 
TBI 1,1) =TBX( 1, J) 

Y = FLOAT( J-1) ’^DY 

IF! ABS (X-XK J) ).LT.0.00001.AND.ABS(Y-Y1I I ) ) .LT. 0.00001) T{ I, J)=TB( 
II ,1) 

IF( ABS( X-XII J) ) .LT.C. 00001. AND. A8S (Y-Y2( I ) ) .LT. 0,00001) T( I ,J )=T8( 
1 I tl) 

IF (ABS(X-X2( J) ) .LT. 0.0000 1. AND. ABS (Y-Y2( I ) ) . LT .0 .00001 ) T(I ,J)=TBI 

11 , 2 ) 

495 CONTINUE 


i 


I 


COMPLETION OF AN ITERATION-CHECKS FOR REQUIRED CONVERGENCE 


IFdTER.SQ.l) WRITE(6,509) 

509 FORMAT(lHl) 

WRITE(6,510) ITER, SUM 

510 CQRMATI5X,* ITERATION = * , I 3 , 3X , ' 6RRCR= • , 
1 E12.5) 
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IFt ITER/NTE*NTE.EQ. ITER) GO TO 520 
I F( SUV.LE.SUMM) GO TO 520 
IF( ITER.GE. I^<AX) GO TO 600 
ITEP.= ITER+1 
GO TO 1 

520 CONTINUE 
WRITE(6,540) 

540 F0RMAT(1HI//////9X, 'TEMPERATURE DISTRIBUTION IN A HOLLOW BLADE • 
I'////// ) 

WRITE (6,541 J DXtDYtOME, XK , NX , NY , SUMM , N 1 , N2 , N3 ,N4,NX0,NXI 

541 F0RMAT{3X,'DX = ' , F S .5 , 5X , ' DY = * , F 8 .5 , 5X , • 0NEGA= ' , F 8 .5 , 5 X , ' THRM .CON 
1D. = *,F8.5/3X,'NX =',I8,5X,'NY =• , I 8 , 5X , • MA X . ERROR* • , ^8 . 5 /3X , * N1 =• 
2,ia,5X,*N2 *',I8,5X,'N3 =',I8,5X,'N4 = ' , I 8/ 3 X , ' NXO * • , I 8 , 5 X, • NX I = 
3* , I6/1H1) 

WRITE! 6,522) 

522 FORMAT!//, 10X,'T-MATRIX, GIVES ALL INTERIOR PO I NTS . I N I T I AL SETTING 

1 OF T IS 1600.0'/10X,' ONLY THOSE INSIDE BLADE ARE CHANGED AND WILL 

2 BE DIFFERENT FROM 1600.00') 

WRITE(6,523) 

523 FORMAT!/, 2X,'J =' , 4X , ' 1 ' , 7 X, ' 2* , 7X , ' 3 ' , 7X , ' 4 ' , 7X , • 5 ' , 7X , ' 6 • , 7X , '7 
1' ,7X,'8',7X,'9',6X,'10',6X,' ll',6X,' 12',6X- '13',6X,'14' ,6X,'15'/ 
22X,' I* ) 

IF(NY.GT.15) WRITE!6,521} 

521 FORMAT! 3X,' 16' ,6X , ' 17' ,6X , ' 1 8* ,6X, ' 19» ,6X, ' 20 ' , 6X , • 2 1 ' , 6X , • 2 2 ' , 6X 
1,*23',6X,'24',5X,'25') 

DO 603 1=1, NX 

603 WRITE!6,525) I , ! T ! I , J ) , J = 1 ,N Y ) 

525 F0RMAT!/,2X, I2,1X, 15F8 , 2 ,2 ! / 5X , 1 5 F 8. 2 ) ) 

WRITE! 6, 524) 

524 FORMAT! 'I',//, 5X, 'TEMPERATURE DISTRIBUTION ON THE BOUNDARY',////) 
WRITE!6,530) NXl 

530 FORMAT!//, lOX, 'BOUNDARY POINTS ON I-LINES =CR LOWER SURFACE OF 

1 OUTER B0UNDRY»/,12X,' (START FROM 1= 2 TO I=',I3,' )') 

WRITE!6,525) ( T8 ! I , 1 ) , I = 1 , NX 1 ) 

WRITE(6,53l) NXl 

531 FORMAT! //, lOX, 'BOUNDARY POINTS ON I LINES FQR UPPER SURFACE OF 

1 OUTER aCUNDPY'/, 12X, •( start FROM 1= 2 '^0 i=*,I3,' }•) 

WR1TE!6,525) ! TB ! I , 2 ) , I = 1 , NXl ) 

IF!N1.EC.N2.AND.N3.E0.N4) GO TO 87 
WRITE (6,532) N11,N21 

532 FORMAT !//,10X, 'BOUNDARY POINTS ON I LINES FOR LOWER SURFACE OF 

1 INNER BOUNDARY'/, 12X, •( START FROM I=',I3,' TO I = ',I3,' )') 

WRITE(6,526) !TB! 1,3 ),I=N11,N21) 

WRITE!6,533) Nll,N21 

533 format !//,10X, 'BOUNDARY POINTS ON I LINES FOR UPPER SURFACE OF 

1 INNER BOUNDARY'/, 12X, ' I START «=ROM I = ',I3,' TO I=',I3,» )•) 

WRITE! 6,526) (TR! 1,4) , I=N 11, M2 1) 

87 CONTINUE 

WRITE (6,534) NYl 

534 'FORMAT!//, lOX, 'BOUNDARY POINTS ON J LINES FOR NEARER SURFACE OF 
I OUTER BOUNDARY' /,12X, '( STAOT FROM J= 2 TO J=',I3,' )•) 

WRITE!6,526)!TBX!1,J),J=2,NY1) 

WRITE(6,535) NYl 
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535 FORMAT (// ,10X, '3CUNDARY POINTS ON J LINES FOR FARTHER SURFACE OF 
1 OUTER BOUNDARY'/, 12X, ' (START FROM J= 2 TO J=',I3,* >•) 

WRITE (6,526) ( T3X( 2 , J ) , J= 2 , NY 1 ) 

IF(N1.EC.N2.AND.N3.EQ.N4) GO TO 89 
WRITE! 6,536) M31,N41 

536 FORMAT!//, lOX, 'BOUNDARY POINTS ON J LINES FOR NEARER SURFACE OF 
IINNEP BOUNDARY ' / 1 2X , * ( ST ART FROM J=',n,' TO J=',I3,' )•) 

WRITE! 6,526) (TBX( 3, J) ,J = N31,N41) 

WRITE(6,537) N31,N41 

537 FORMAT!//, lOX, 'BOUNDARY POINTS ON J LINES FOP FARTHER SURFACE OF 
lINNER BOUNDARY '/ 12X ,»( START FROM J=',I3,' TC J=',I3,» )•) 

WRITE! 6,526 ) !T 8X ! 4, J ) , J= N3 1 , N4 1 ) 

WRITE!6,401) 

401 FORMAT !1H1) 

526 F0RMAT(3!/,5X,15F8.2)) 

89 CONTINUE 

WRITE !6,801) 

DO 701 J1=2,NY1 
DO 702 1 = 1, NX 
J=NY+l-Jl 

XR=FLOAT! I-l )MOX j 

IF! !XR-X1! J) ) .GE.O, 0001 ) 1 1 = 1 i 

IF(ABS(XR-X1!J)).LE. 0.0001) 11=1+1 j 

IF! !XP-X1!J)). GE.O, 0001) GO TO 704 

702 CONTINUE ! 

704 DO 703 1=11, NX 
XR=FLOAT! I-l )*0X 

IF! (XR-X2U )) .GE.O. 0001 ) 12=1-1 
IF!ABS!XR-X2!J) ) .GE .! 0 .9999*0 X ) ) 12=12-1 
IF! (XR-X2!J)). GE.O. 0001) GO TO 705 

703 CONTINUE 

705 CONTINUE | 

J=NY+1-J1 j 

IF!J.E0,2 / WRITE(6,824) J ,TBX! 1,J ) , ( T! I, J) , 1 = 11, 12) ,TBX!2, J) j 

IFIJ.EQ.3 ) WRITE!6,823) J,TBX! 1, J ) ,(T( I, J) , 1=11, 12) ,T8X!2, J ) 

IF!J.E0,4 ) WRITE! 6,322) J ,TBX! 1 ,J) , !T( I ,J ), 1=11,12) ,TBX!2,J) 

IF!J.E0.5 ) WRITE(6,821) J, TBX! I ,J ) ,( T{ I, J) , 1=11, 12) ,TBX!2, J) 

TF(J,EQ.6 ) WRITE!6,820) J,T8X!1,J) ,!T! I,J) , 1=11, 12) ,TBX!2,J) 

IF!J,E0.7 ) WBITE!6, 819) J ,TBX ! 1, J) , (T( I , J) , 1 = 11, 12) ,TBX!2, J) 

IF!J,E0.8 ) WRITE !6 ,818) J ,TRX ( 1 , J ) , !T! I , J ) , 1 = 11, 12) ,TBX!2, J) 

IP!J,EQ.9 ) WPI TE!6,817) J ,TBX( 1 , J) , !T! I , J) , 1 = 11, 12) ,T8X!2, J) ■ 

IF(J.EO.IO) WRITE!6,816)J,TBX!1,J),(T!I,J) , 1=11,12) ,T3X!2,J) j 

IF (J.EQ.ll) WRITE (6 ,815) J, TBX! 1, J) , !T( I , J ) , I = 1 1 , I 2 ) , T 3X ( 2 , J ) 

IF(J.EQ.12) WRITE!6,ai4) J,TBX( l,Jj ,!T( I,J) ,1=11,12) ,T8X!2, J) j 

IF(J.EQ.13) WRITE!6,81 3) J,TOX( 1, J } , (T( I , J) , 1=11, 12) ,TBX(2, J ) i 

IF(J,EQ.14) WRITE!6,812) J, TBX! 1 , J) , !T ! I , J) , 1=11, 12 ) .TBX (2, J) 

IF!J.EQ.15) WRITE!6, 811) J,T3X! 1,J), IT! I, J) , 1=11,12) ,TRX(2,J) | 

IF!J.E0.16) WRITE !6,ai0) J, TBX! 1 , J ) , !T( I , J ) , 1=11, 12) ,TBX!2, J ) 

IP! J. £0,17) WPITE(6,809) Jt TBX! l,J),!T!I,Jj ,I=I1,I2),TRX(2,J) 

IF!J.E0.18) WRrTE(6,30 8)J,TBX(l,J),(T(I,J), 1 = 11,12) ,TBX(2, J) 

IF ! J.E0.19) WRITE !6 ,807) J, TBX! 1, J ) , !T( I , J) , I = IL, 12) .TBX(2, J) 

IF!J.EQ.20) WRITE!6,806)J,T3X!1,J),!T! ItJ) ,I=Il,IZ: ,T8X!2,J) 

IFU.EQ.21) WPITE(6,80 5)J,T3X!1,J),!T(I,J),I=I1,I2),TBX(2,J) 



j 

} 
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IF( J.EQ.22) WRITE(6,804) J,TBX( ltJ)f (T(I, J) , 1=11, 12) .TBX(2t J) 
IFU.EQ.23 ) WRITE(6,e03) J,TBX( 1 , J ) , ( T ( I , J ) , I = 1 1 ♦ 12 ) , TBX ( 2 , J ) 
IFU.EQ.24) WO ITE 16,80 2) J , TBX( lrJ),(T(I,J),I = Il,I2) , TBXI2, J) 

801 FORMATl 1H1///,15X, 'TEMPERATURE LxSTRIRUTION INSIDE THE BLADE AND 0 
XN ITS SURFACE'///2X, 'J=* ) 

002 F0RMAT(/2X,I2, 4X,20F8.21 

803 FORMAT{/2X,I2, 6X,20FB.2) 

804 F0RMAr(/2X,I2, 8X,20F3.2) 

805 FORMAT (/2X, 12, lOX, 13F8 .2/24X, 7F3.2 ) 

806 F0RMAT(/2X,I2,12X, 13F8 .2/24X , 7F3. 2 ) 

807 F0RMATI/2X, 12, 14X , 1 3 F8 .2 /24X , 7F8 . 2 ) 

808 FORMAT! /2X,I2,16X, 13F8.2/24X, 7F8.2 ) 

809 FORMAT ( / 2X , I 2 , 18X, 1 3F8 . 2 /24X , 7^8 . 2 ) 

810 FORMAT(/2X,I2,20X,20F8.2) 

811 FORMAT! /2X, I2,22X,20F8.2 ) 

812 F0RMAT!/2X,I2,24X,20F3.2) 

813 FORMAT!/2X,I2,26X,20'=8.2) 

814 FQRMAT!/2X,I2,28X,20F8.2) 

815 PQRMAT!/2X,I2,30X,20F8.2) 

816 FORMAT!/2X,I2,32X,20F8.2) 

817 =QRMAT! /2X, I2,34X,20F8.2) 

818 FORMAT !/2X, I2,36X,20F8.2 ) 

819 FORMAT! /2X,I2,38X,20F8.2) 

820 F0RMAT!/2X,I2,40X,20F8.2) 

821 FORMAT!/2X,I2,42X,20F8.2 ) 

822 FORMAT! /2X, I2,44X,20F8.2) 

823 FORMAT(/2X,I2,46X,20F8.2) 

824 FORMAT! /2X, I2,48X,20F8.2) 

701 CONTINUE 

IFISUM.LE.SUMM) GO TC 600 
I TER = ITER+1 
GO TO 1 
600 CONTINUE 
633 WRITE!6,665) 

665 FORMAT! 1H1,38X,' ISOTHERMAL LINE LOCATIONS * / ,9X, ' I ' ,4X, * J* ,3X , • T • , 
14X,'T - 1*,6X,'T - 2 • ,8X,'FRAC»,/) 

DO 601 1=2, NXl 
IA=IFIX!Y1! I)/DY+0.0001) + 1 
I8 = IFIX!Y3! I)/OY+0.000 1)+2 
IC=IFIX!Y4!I )/0Y+0.0001) +1 
I0=IFIX!Y2! I )/DY +0.0001) 1-2 
I8=IB-1 
ID=I0-1 

IF! I.GT.M.AND.I.LT.NZ) GO TO 629 
IT1=T0! 1,1) 

IT2=TB! 1,2) 

IL=IT1 

IH=IT2 

IF! IT1.GE.IT2) IL=IT2 
IF ! ITl.GE, IT2) IH = ITl 
I AA = IA 

DO 623 J=IAA,ID 

IF! Tt I , J ).LT. FLOAT! ID ) I L = I F I X ! T ! I , J ) ) 
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IF(T(I,J).GT.FLOAT( IH) )IH=I=IX{T( I,J)) 

620 CONTINUE 

623 CONTINUE 

DO 624 II=IL,IH 
DO 625 J=IA,ID 

IFHFLOATIJ )#DY) .GT.Y2I I n T ( I , J 4- 1 ) = T3 ( 1 , 2 ) 

IF { Yin !.GT. {FLOAT! J-1 )=!<0Y) ) T ( I ♦ J ) = TB (I ♦ 1 ) 

TV = FLOAT{ 1 1 ) 

IF(T( I , JI.LT.TV.ANO.Td, J+D.GE.TV) GO TO 627 
IF(T(I,J).GE.TV.AND.T( I,J + 1) .LT.TV) GO TO 627 
GO TO 625 
627 CONTINUE 

RX = {T{ I ,J >-TV)/( TCI , J)-T( I tJ + 1) ) 

WRITE [6, 666) I,J»II,T{I,J),T(I»J+1),RX 
625 continue 

624 CONTINUE 
GO TO 601 

62*3 CONTINUE 

IT1=TB(I ,1) 

IT2=TB( 1,3) 

IT3=TB( 1,4) 

IT4=TB( 1,2) 

IL=IT1 

IH-IT2 

ILl=IT3 

IHL=IT4 

IF(IT1.GE.IT2) IL=IT2 
IF{ IT1.GE.IT2) IH=IT1 
IF{ IT3.GE. IT4) ILl = IT4 
Icnr3.GE.IT4) IH1=IT3 
DO 604 n=IL,IH 
DO 605 J=IA,IB 
TV = FLOAT( 1 1 ) 

IFIYK I ).GT. (FLOAT ( J-1 )#0Y) ) T ( I , J ) = TB ( I , 1 ) 
IF((FLOAT(J ) #DY) .GT.Y3 ( I ) ) T ( I , J+-1 ) =TB ( I , 3 ) 
IF(Tn,J).LT.TV,AND.T( I,J+1) .GE.TV) GO TO 607 
Ic(Tn,J).GE.TV.AND.T( I,J + 1).LT.TV) GO TO 607 
GO TO 605 

607 RX=(T( I,J )-TV)/(T(I ,J)-T( I,J<-1) ) 

WRITE! 6,666) I , J, I I ,T( I , J) ,T!I ,J-H1 ,RX 
666 F0RMAT(5X,3I5,3F11.4) 

605 CONTINUE 
604 CONTINUE 

00 614 II = IL1,IH1 
DO 615 J=IC,ID 
TV-FLOAT! I I ) 

I=!Y4!I ) .GT. (FLOAT! J-1 )*DY) ) T ! I , J ) = T9 ! I , 4 ) 

IF ( (FLOAT! J-1) *OY).GT. Y2( I ) ) T(I , J + 1 ) =TB ( I , 2 ) 

1 F(T( I , J) .LT.TV.ANO.T! I , j+l) .GE.TV) GO TO 617 
IF(T!I,J).GE,TV.ANO.T! I,J+1),LT.TV) GO TO 617 
GO TO 615 

617 RX=!T! I,J)-TV)/(T!I,J)-T(I,J+1)) 

WR ITE( 6,6 66) I , J , II , T ( I , J ) , T ( I , J + 1 ) , R X 
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615 CONTINUE 
614 CONTINUE 
601 CONTINUE 
645 CONTINUE 
STOP 
END 
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SUBROUTINE M ESH ( I , J » TNEW ) 

COMMON T(40,25) ,TB(40»4) ,TBX(4,25) tTG(40,4) ,TGX(4,25),H(40,4),HX(4 
1, 25) ,Y1(40) , Y2(40 ) , Y3( 40) ,Y4(40) ,Xlt 25) ,X2( 25) ,X3( 25) »X4(25) ,BX( 12 
20) tBY( 120) ,BXI (120) fPYK 120) , IP( 120) 

COMMON DX,DY,NX,NY,NXO,NXI ,ITER»NltN2»N3»N4 


SUBROUTINE FINDS S I- 1 , S I -2 t OEL-U , DEL-2. FOR EACH INTERIOR POINT AND 
GETS NEW TEMPS FOR EVERY INTERIOR POINT 


PAR=FLQAT( I-1)*DX 

IF( (PAR-XK J)) ,LE.DX)S1= (PAR-Xl ( J ) )/DX 
IFt ( PAR-Xl U ) ) .LE,DX)T1=TBX( 1,J) 

IP ( (PAR-XK J) l.LE.CX) GO TO 20 
IF ( J,GT.N3.AND,J.LT.N4)GC TO 10 
GO TO 15 

10 IF (PAR,GE.X4(J) )G0 TC 12 
GO TO 15 

12 IF( (PAR-X4U) ).LE.DX)S1=(PAR-X4( J) )/DX 

IF( (PAR-X4( J )) .LE.CX) T1=T0X(4,J) 

IF KPAR-X4( J) ).LE.DX) GO TO 20 
15 Sl=1.0 

Tl=T( I-l, J) 

IF(A3S(X1 (J)-PAR+DX) .LT. 0.00001) T1=TBX( I , J) 
IF(ABS(X4( J)-PAR + DX ) ,LT, 0.00001) T1=TBX( 4,J ) 
20 CONTINUE 

IF( (X2( J)-PAR) .LE.DX)S2=IX2( J)-PAR )/DX 
IF ( (X2(J)-PAR) .LE.DX ) T3=TBX(2,J ) 
TF({X2(J)-PAR) .LE.DX) GO TC 35 
IF( J.GT.N3.AND.J.LT.N4)G0 T025 
GO TO 30 

25 IF(X3(J).GT.PAR)GC TC 27 

GO TO 30 

27 IF(( X3( J)-PAR),LE.DX)S2=(X3(J)-PAR)/DX 

IF((X3(J)-PAR),LE.OX)T3=TBX(3,J) 

IF((X3( J)-PAR). LE.DX) GO TO 35 
30 S2=l. 

T3 = T( I + l, J ) 

IF (ABS(X2 ( J)-PAR-DX ) .LT. 0.00001 ) T3=T8X( 2» J ) 
IF (ABS(X3( J) -PAR -OX). LT.O. 00001) T3=TBX(3,J) 
35 PAR=FLOAT( J-l)*OY 

IF{ (PAR-YK I)),LE.0Y)D1=(PAR-Y1{I))/DY 
I<= ( (PAR-YK I ) ) .LE.DY) T2=Tfi( 1,1) 

IF( (PAR-Yl { I ) ) .LE.DY) GO TO 50 
IF{ I.GT.M.AN0.I.LT.N2)G0 TC40 
GO TO 45 

40 IF(PAR.GE.Y4( I) )G0 T042 
GO TO 45 

42 IF(( PAR-Y4( I ) ) .LE.CY) D1 = (PAR-Y4( I ) ) /DY 
IF((PAR-Y4( I)),LE.CY)T2 = TR(I,4) 

IF { (PAR-Y4I I ) ' LE.DY) GO TC 50 
45 Dl=1.0 
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T2=T{ I ,J-1) 

IF (ABS( YU I )-PAR+DY ).LT. 0.0000 1 ) T2 = TB( I ,1 ) 
IP(ABS(Y4( D-PAR4-DY) ,LT, 0.00001) T2=TB(I i4) 
50 CONTINUE 

IF ( (Y2(I )-PAR) .L£.DY)D2=(Y2 ( I )-PAR) /DY 
IF{{Y2(1 )-PAP, ) .LE.DY) T4=T8( It2) 

IF( (Y2( I)-PAR) .LE.OY) GO TO 65 
IF( I.GT.Nl.AND. I.LT,N2)G0 TC55 
GO TO 60 

55 IFIY3( I) .GT.PARIGO T057 
GO TO 60 

57 IFMY31 I )-PAR ) ,L E . CY ) 02= ( Y3 (I )-P AR ) /DY 
IFnY3II )-PAR).LE.DY)T4 = TB{ 1,3) 

IF( (Y3( I )-PAR) .LE.DYJ GO TO 65 
60 02=1.0 

T4=T( I ,J+1) 

IF (ABS(Y2( I )-PAR-DY) .LT.0.00001) T4=TB( 1,2) 
IF (ABS{Y3 (I )-PAR-DY) .LT. 0.00001) T4=TB( I ,3) 
65 CONTINUE 

A1=T1/S1/ (S1 + S2 ) 

A2=T2/D1/ID1+D2) 

A3=T3/S2/{ S1+S2) 

A4=T4/D2/ (D1+D2) 

E =1. /SI/ S2+{DX/DY)*(0X/DY)/ 01/02 
TNEW= {A1 + A3 + DX=0X/(0Y>!'0Y) =«'( A2 + A4) )/E 
RETURN 
END 
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SUBROUTINE S LO ( X I , Y 1 , X2 , Y 2 , X3 , Y3 , A ♦ B , C ) 

INOS EQUATION THRU THREE POINTS AS Y= A+B$ X+C^X«X 


CALL UNDFLW 

IF (ABS( Xl-X2).LE.O.Obl.ANn,ABStX2-X3 I.LE.O.OOl ) GO TO 10 

XS3=X3*X3 

XS2=X2*X2 

XSl=Xl*Xl 

Dl=X2*XS3-X3*XS2 

D2=X1*XS3-X3^XS1 

03=X1»XS2-X2*XSI 

0=01-02+03 

IF(A3S(D).LE. 0.000001) GO TO 10 
Al=Yl*(X2*XS3-X3*XS2) 

A2=Y2=*«( X1*XS3-X3*XS1 ) 

A3=Y3*( X1*XS2-X2^XS1 ) 

A=( A1-A2+A3) /O 
ei=Y2’!'XS3-Y3*XS2 
B2=Y1«XS3-Y3*XS1 
B3=Yl=i'XS2-Y2*XSl 
B=( B1-B2+B3) /O 
C1 = X2*Y3-X3=*=Y2 
C2=X1*Y3-X3*Y1 
C3=X1*Y2-X2*Y1 
C=(Cl-C2+C3)/0 
RETURN 
10 A=100. 

B=100. , 

C=100. 

RETURN 

END 
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SUBROUTINE CUR(B2tC29BtC, I 1 

COMMON T(40t25) ,TB ,TBX(4,25) ,TG(40,4) ,TGX(4f25),H(40,4)f HX{4 

1»25) ,Y1 (40) ♦ Y2140) ,731 40) ,Y4(40» ,X1( 25) »X2( 25) ,X3l 25) ,X4(25) ,BX( 12 
I 20) ,BY(120) ♦BXH120) ,BYI ( 120) » IP ( 120) 

! COMMON DX,DY,NX,NY,NXC,NXItITER,Nl,N2,N3,N4 


subroutine to find adjacent two POINTS FOR THE PQINT(82»C2) AND THEN 
PASSES A CURVE Y = A+BX<'CX2 AND DETERMINES A,B,C. B,C COEFFICIENTS 
ARE RETURNED TO MAIN PROGRAM TO GET SLOPE OF CURVED 8CUN0ARY AT THE 
POINT 

. i ' 


IF( I .EQ.l )G0 TO 25 
DO 10 M^ltNXO 

IF(ABS(B2-BX{M)). LE, 0,0001. AND, A6S{C2-BY(M) ).LE. 0.0001) GO TO 12 
10 CONTINUE 

12 IF (M.EO.l ) GO TO 15 \ 

Bi=BXlM-l) I 

C1=BY(M-1) j 

GO TO 13 

15 B1=8X(NX0) I 

C1=8Y(NXC) I 

13 CONTINUE ; 

83 = BX(M+1) 

C3 = 3YIM+1) 

IFIM.FQ.NXO) 83=8X11) 

IF(M.EO.NXO) C3=BY(1) 

IF(B3.EQ.B2-AND,C3.E0.C2)GG TO 14 
GO TO 37 
14 83 = BX(M+2) 

C3 = BY(M+2) 

GO TO 37 

25 no 30 M = l,NXI 

IF(B2.EQ.0XI (M),AND.C2,EQ.BYI[M) )GG TO 32 
30 CONTINUE 

32 IF(M.EC.l) GO TO 34 
BI=BXI(M-1) 

CI=BYIIM-1) 

IF(M.E0.17) Bl=5.8 
IF(M.EQ.17) Cl=1.74 
GO TO 33 

34 31=BXItNXl) 

C1=BYI ( NXI) 

33 CONTINUE 
B3=BXI(M+1) 

C3=BYI(M+1) 

ip(M.6Q.Nxi) B3=exim i 

IF(M.EO.NXI) C3 = 8Yim i 

I=( F^.EQ.02.AND.C3.EO.C2)GO TO 35 i 

GO TO 37 
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1 


I 


IV G LEVEL 21 CUR DATE = 743^6 

35 B3=BXI { f^+2) 

C3=BYI (M+2) 

37 CALL SLO(Bl,Cl,B2,C2,33,C3,A,e,C) 

IP{ ITER.EO.l) WRITE(6t40n,M,Bl,Cl,B2tC2fB3,C3,A,B,C 
40 FORMAK 3 X , 12 » I 5 , 9 F 1 0 .5 ) 

RETURN 

ENO 
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